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Abstract: 

Charged asymptotically AdSc, black branes are sometimes unstable to the condensation of charged 
scalar fields. For fields of infinite charge and squared mass —4 Herzog was able to analytically determine 
the phase transition temperature and compute the cndpoint of this instability in the neighborhood of 
the phase transition. We generalize Herzog's construction by perturbing away from infinite charge in 
an expansion in inverse charge and use the solutions so obtained as input for the fluid gravity map. 
Our tube wise construction of patched up locally hairy black brane solutions yields a one to one map 
from the space of solutions of superfluid dynamics to the long wavelength solutions of the Einstein 
Maxwell system. We obtain explicit expressions for the metric, gauge field and scalar field dual to an 
arbitrary superfluid flow at first order in the derivative expansion. Our construction allows us to read 
off the the leading dissipative corrections to the perfect superfluid stress tensor, current and Josephson 
equations. A general framework for dissipative superfluid dynamics was worked out by Landau and 
Lifshitz for zero superfluid velocity and generalized to nonzero fluid velocity by Clark and Putterman. 
Our gravitational results do not fit into the 13 parameter Clark-Putterman framework. Purely within 
fluid dynamics we present a consistent new generalization of Clark and Putterman's equations to a set 
of superfluid equations parameterized by 14 dissipative parameters. The results of our gravitational 
calculation fit perfectly into this enlarged framework. In particular we compute all the dissipative 
constants for the gravitational superfluid. 
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1. Introduction and Summary 

It was pointed out by Gubser that charged asymptotically AdS^ black branes are sometimes unstable 
in the presence of charged scalar fields. The endpoint of this instability is a hairy black brane: a black 
brane immersed in a charged scalar condensate. The AdS/CFT correspondence maps the hairy black 
brane to a phase in which a global C/(l) symmetry is spontaneously broken by the vacuum expectation 
value of a charged scalar operator. In condensed matter physics a phase with a spontaneously broken 
global U{1) symmetry is referred to as a superfiuid. 

In this paper we study aspects of the fiuid dynamical description of superfluids using the AdS/CFT 
correspondence. Our work builds on a large body of earlier studies (see e.g. ^, |[ |[ ^, 0, ||, ||, ^ 
[l], [I 0, 0, [8[ 1^, H) but differs from these works in its emphasis on the study 

of dissipative terms in super fluid dynamics. 

As we describe in detail in the next section, the variables of relativistic superfluid dynamics consist 
of two velocity flelds; the normal fluid velocity u'' and a superfluid velocity fleld m^; together with a 
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temperature and chemical potential field. The superfiuid velocity is the unit vector in the direction of 
—£,fj. where is the gradient of the phase of the scalar condensate. Conservation of the stress tensor 
and charge current together with the assertion that is curl free constitute the equations of superfiuid 
dynamics. These equations constitute a closed dynamical system once they are supplemented with 
constitutive relations that express the the stress tensor, charge current and the component of along 
the normal velocity as functions of the fluid dynamical variables. As superfiuid dynamics is a long 
distance effective field theory, it is natural to specify the relevant constitutive relations in an expansion 
in derivatives. 

Over fifty years ago Landau and Tiza [|3[ presented a simple and elegant proposal for the struc- 
ture of superfiuid constitutive relations at leading (zero) order in the derivative expansion. Landau 
and Tiza (see below) proposed a form for the constitutive relations that is entirely determined by a 
single thermodynamical 'Free Energy' P(T, fi, ^). An obvious first question in the study of holographic 
superfluids is the following: do the perfect fiuid constitutive relations of holographic superfiuids re- 
spect the Landau-Tiza ansatz? This question was largely answered in the affirmative in a beautiful 
recent paper by Sonner and Withers 0. These authors used Einstein's equations and the holographic 
dictionary to demonstrate that the stress tensor and charge current of a homogeneous stationary holo- 
graphic superfiuid fiow takes the form predicted by Landau and Tiza, with the free energy or pressure 
interpreted as the value of Einstein's action for the relevant bulk solutions. We re-derive and slightly 
extend the results of Sonner and Withers in §^ and Appendix ^ below. As we review in Appendix 
^ a beautiful feature of the gravitational derivation of the Landau-Tiza model is its abstract nature. 
The results of Sonner and Withers are derived on general grounds, and do not use the explicit form 
of the gravitational solution dual to homogeneous stationary superfiuid flows. This is fortunate as 
no completely explicit analytic solutions are known for static hairy black branes ( i.e. holographic 
superfluids at rest) much less for hairy black branes in motion. 

There also exists a large literature on the subject of dissipative corrections to the equations of 
superfluid dynamics, accurate to first order in the derivative expansion (see e.g. |25[ ^ |2^, |o[ 
^ ^ ^ |4[ m |6[ ^ ^). In particular, a general first order theory of dissipative fiuid dynamics 
was presented by Landau and Lifshitz (under the assumption of small superfiuid velocities) and 
was generalized to apply to flows with finite superfluid velocities in |6| . It turns out that the most 
general one derivative corrections to the equations of perfect superfluid dynamics are parameterized 
by 36 dissipative corrections (assuming parity invariance). Like Landau and Lifshitz, in [ ps] , the 
authors assumed that the entropy current of superfluid dynamics takes a natural 'canonical' form at 
one derivative order. Once again, like Landau and Lifshitz, the authors in |25| p6| used the requirement 
of local positivity of entropy production, together with the 'Onsager Reciprocity relations' to cut the 
parameter space for physically permissable equations of superfluids down to a 13 parameter set (see 
Appendix VI in ) . To our knowledge this Clark-Putterman formalism has not been tested 

by any first principles dynamical calculations within a quantum field theory (e.g. has not been derived 
from a field theoretically motivated set of Boltzmann equations). The AdS/CFT correspondence offers 
us the opportunity to test this abstract formulation of viscous superfluid dynamics. 

In this paper we perform such a test in a very particular holographic super fluid. The system we 
study is a small perturbation of the 'analytic superconductor' studied by Herzog in In very broad 
terms we use the fluid gravity map |39[ Q ^ to derive the equations of superfluid 

dynamics from Einstein equations to first order in the derivative expansion. We then read off the 
constitutive relations of the stress tensor, charge current and Josephson equation at first order in the 
derivative expansion. In the rest of this introduction we describe our calculations, their results and 
their interpretation in more detail. 
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As we have mentioned we work with a generahzation of the model studied by Herzog Herzog's 
model consists of the Einstein Maxwell system interacting with a charged scalar field of = —4 in 
the so called probe limit of infinite charge e. Herzog demonstrated that this system undergoes a second 
order phase transition towards superfiuidity whenever |^| > 2. The stable gravitational solution, for 
1^1 just larger that 2, has a background scalar vev. Let e denote the value of this vev. In Herzog 
analytically determined the relevant bulk solutions perturbatively in e and separately in the difference 
between superfluid and normal velocities. 

To start with we generalize Herzog's infinite charge solutions beyond the strict probe approxima- 
tion, to first nontrivial order in the ^ . This generalization is necessary in order to allow for the study 
of the response of the normal velocity and temperature fields to the dynamics of the superfiuid velocity 
and chemical potential fields. As a check on our algebra we explicitly verify that our solutions obey 
all the predictions of the Landau Tiza model, to the order that we are able to compute, in accordance 
with the results of Sonner and Withers and of Appendix We then proceed to use these solutions 
as raw ingredients for the fluid gravity correspondence |^ ^ p9| , |4^ , ^ . 

Following the procedure of the fluid gravity correspondence, we search for solutions of the Einstein 
Maxwell scalar system that tube wise approximate the stationary solutions described in the previous 
paragraph. More explicitly, we study a perturbative expansion to the solutions of Einstein's equations 
whose first term is given by the stationary solutions of the previous paragraph -written in ingoing 
Eddington Finklestein coordinates - with the eight parameters of equilibrium superfiuid flows replaced 
by slowly varying functions of spacetime. The configuration described in this paragraph does not obey 
the bulk equations; however it may sometimes be systematically corrected, order by order in boundary 
derivatives, to yield a solution to these equations. This procedure works if and only if our eight fields 
are chosen such that £,'^{x) is curl free, and such that the energy momentum and charge current built 
out of these fields is conserved. The constitutive relations that allow us to express the stress tensor 
and charge current in terms of fluid dynamical fields is generated by the perturbative procedure itself. 
In other words the output of our perturbative procedure is a set of gravitational solutions that are 
in one to one correspondence with the solutions of superfiuid dynamics, with superfluid constitutive 
relations that are determined by the bulk gravitational equations. 

Note that the construction described in the previous paragraph is carried out in a triple expansion. 
We follow Herzog to expand our equilibrium solutions in a power series in the deviations from criticality 
(let us denote the relevant parameter by e) , and further expand these solutions in a power series in ^ . 
We then go on to use the solutions as ingredients in a spacetime derivative expansion. We would like 
to emphasize that this procedure is sensible only if the derivatives times mean free path are assumed 
to be parametrically smaller than e. The physical reason for this is that precisely at e = we have 
a new massless mode, corresponding to fluctuations of the vev of the charged scalar field. When e is 
nonzero this mode is no longer massless, but it is light at small e. The fluid dynamical description 
ignores the dynamics of this light mode. This is justified only when derivatives are all much smaller 
than the mass of this mode, i.e. when derivatives are parametrically small compared to e. The fluid 
dynamical expansion breaks down if derivatives are held flxed as e is taken to zero. This fact formally 
shows up in the blow up of several dissipative fluid coefficients at small e in the equations presented 
in this paper. 

We have implemented the perturbative procedure that determines the gravitational solution dual 
to superfluid dynamics to flrst order in the derivative expansion, and thereby determined the stress 
tensor, charge current and Josephson equation of our holographic superfluid to first order in the 
derivative expansion. Our results pass several consistency checks and also have several attractive 
features. To start with our results respect Weyl invariance. This fact, while necessarily true of the 
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results of any gravitational calculation in an asymptotically AdS space, is not algebraically automatic, 
and so yields consistency check on the rather involved algebraic manipulations of our paper. 

More importantly, we are also able to compute the natural gravitational entropy current for our 
gravitational fluid flow. As explained in , a gravitational construction of fluid dynamics is always 
automatically accompanied a local entropy current that is guaranteed, by the area increase theorems 
of general relativity, to be of positive divergence. We demonstrate by explicit computation in our 
particular solution that the gravitational entropy current of p9| agrees precisely with the 'canonical' 



entropy current (see §2.5 for details) of Clark and Putterman (and Landau Lifshitz) [g5|, |27|. We 
regard this agreement as a nontrivial check of one of the central physical assumptions of the Landau 
Lifshitz and Clark Putterman formulations of dissipative fluid dynamics. ^ Despite this fact, however, 
the results we obtain for the dissipative coefficients for the gravitational superfluid described in this 
paper do not fit into Clark-Putterman's 13 parameter framework for dissipative fluid dynamics ^ . 

Given the fact that the gravitational entropy current agrees with the current employed by Clark 
and Putterman [|5[ |6| , it is puzzling that the gravitational results do not fit into Clark and Putter- 
man's framework. We believe that the resolution to this puzzle is that Clark and Putterman missed 
a parameter in their analysis. In §^ below we demonstrate that the most general equations of super 
fluid dynamics consistent with positivity of the canonical entropy current is parameterized by 21 dis- 
sipative coefficients rather than 20 as mistakenly asserted by Clark and Putterman. Imposition of the 
7 Onsager relations then leaves us with a 14 (rather than 13) parameter set of consistent equations of 
first order dissipative super fluid dynamics. We present this generalized 14 parameter set of equations 
for dissipative fluid dynamics in §^ below. It turns out that our gravitational results fit perfectly into 
this enlarged 14 parameter framework In §H below we explicitly list the dissipative parameters for 
our gravitational superfluid. 

Our paper is organized as follows. In §^ and §^ we present a general framework for relativistic 
superfluid dynamics including dissipative effects at first order in the derivative expansion. These 
sections are purely fluid dynamical and makes no reference to the equations of gravity or the AdS/CFT 
correspondence. The chief new result of §^ is the presentation of the 14 parameter set of equations 
that generalize Clark and Putterman's 13 parameter equations We also demonstrate that 

the equations of Weyl invariant superfluid dynamics are parameterized by 10 dissipative parameters, 
and present the most general form for these equations. In below we review and slightly generalize 
Sonner and Wither's gravitational derivation of the Landau Tiza 2 fluid model from gravity. In §^ we 
perturb Herzog's construction of hairy black branes away from the strict large charge or probe limit. 
In ^ we use the results of §^ as an input into the fluid gravity map to generate the gravitational 
solutions dual to superfluid flows. We compute the dissipative part of the stress tensor, charge current 
and entropy current dual to superfluid flows and also verify the Weyl invariance of our results. In s|| we 



^We should, however, emphasize that we have checked this agreement only to leading nontrivial order in the e 
expansion. It would certainly be useful to verify this agreement - and all the other results of our paper - at higher orders 
in the e expansion, but we leave that for future work. 

^In an earlier version of this draft we compared our results to those of Landau and Lifshitz but not to those of Clark 
and Putterman, as we were unaware of their work. We reported that our results did not agree with the predictions 
of Landau and Lifshitz; while this is true, the comparison itself is inappropriate, as Landau and Lifshitz apaprcntly 
intended their analysis to apply only to the limit of zero superfluid velocity, while in this paper we work at arbitrary 
superfluid velocity. In this version of our paper we have instead compared our results to the predictions of Clark and 
Putterman who explicitly work at finite superfluid velocity. We once again find disagreement with their predictions, and 
have motivated us to construct a slight generalization of the Clark-Putterman formalism, with which our gravitational 
results now agree. We thank C. Herzog and A. Yarom for making us aware of Clark and Putterman's work. 

'^Imposin g co nformal invariance (which is relevant for the gravitational calculation) sets 4 out of these 14 parameters 
to zero (see §p.^ for more details). 
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transform our gravitational stress tensor and charge current to a frame convenient for the study of fluid 
dynamics, and demonstrate that our gravitational results are a special case of the 14 parameter set of 
dissipative fluid dynamical equations derived in We also explicitly list the values of all dissipative 
parameters of our gravitational superfluid. In appendix ^ we extend the discussion of the canonical 
entropy current in §2.5 to demonstrate that it is independent of the choice of frame and also compute 
its divergence. In appendix ^ we present an abstract geometrical derivation of the analogue of Gibbs 
Duhem relations for superfluids in thermodynamic equilibrium directly from gravity. In appendix 
^ we show the existence of an upper bound for the superfluid velocity beyond which the superfluid 
phase is unstable. Finally, in appendix ^ we complete the discussion of j |7.5| by manifestly recovering 
complete 5*0(3) rotational invariance in the zero superfluid velocity limit. 

Note Added: Our paper has substantial overlap with which was posted on the ArXiv simul- 
taneously with the first version of this paper. 



2. Review of relativistic superfluid dynamics 
2.1 Relativistic superfluids in equilibrium 

Consider a relativistic quantum field theory with a conserved U{\) charge. In the sequel we denote 
the conserved U{1) current of this theory by J'' and the conserved stress tensor by T^^^ . Consider 
this system at a finite temperature T and finite chemical potential /i. It is conceivable that a charged 
operator O (charged under the U{1) charge described above) develops a nonzero vev over a certain 
range of temperatures and chemical potentials. Whenever this happens our system is said to display 
superfluidity. Superfluidity, then, is associated with the spontaneous breakdown of a global U{1) 
symmetry. One of the most striking facts about a superfluid is that it admits more stationary homo- 
geneous solutions to the equations of motion, on i?"^'^, that might naively have been supposed. Any 
system with a U{1) charge admits at least a 5 parameter set of homogeneous solutions on R^'^\ these 
solutions represent the system in equilibrium at temperature T and chemical potential /u, moving at 
a uniform four velocity . However superfluids actually admit an 8 parameter set of homogeneous 
stationary solutions, as we will now explain. 

Let -0 denote the phase of the condensed operator O. ■(/' is effectively a massless scalar fleld in 
the superfluid phase. A massless scalar fleld admits solutions of the form il){x) = e^^x'^ for arbitrary 
constant values of (here e is the charge of the operator O). It turns out that a superfluid admits 
homogeneous stationary solutions at every constant value of T, and f^. 

There is another way to think of the 8 parameter set of solutions listed above. We can 'gauge' the 
global symmetry of the field theory, and so couple the theory to a non dynamical gauge field A^. The 
solutions described above may all, equivalently, be thought of as solutions with constant values of the 
phase ip, but with the non dynamical gauge fields 

where e is the charge of the field The equivalence of these two ways of describing these solutions follows 
from the fact that a phase e^^cc^ of a field of charge e is gauge equivalent to the gauge field listed 
above. 

Finally some definitions. We define 

^l = ^^uf". 

fi is referred to as the chemical potential of the system (or sometimes as the chemical potential of the 
normal part of the system). This terminology is reasonable as /z is equal to the asymptotic value of 
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the time component of the gauge field in the frame in which the normal fluid is at rest (i.e. in the 
frame in which = (1, 0, 0, 0)). Also let 



We define the normalized 4 vector 

and will refer to as the 'superfluid velocity' in the solution. Roughly speaking, is the four 
velocity at which the stuff associated with the condensate, (O) moves. 

With all this terminology in place, we can reword the central assertion of the previous paragraph 
as follows: the stationary non dissipative solutions of a superfluid are labeled by the three components 
of u'^, the 'velocity' of the normal component of the fluid, the three components of it^, the velocity 
of the superfluid (or condensate) of the fluid, the temperature T and the chemical potential /x. Most 
importantly we have exactly stationary solutions in which the 'normal' and 'superfluid' parts of the 
system move at arbitrary speeds with respect to one another. 

We now turn to the quantitative characterization of the stationary homogeneous flows described 
above. We would first like quantitative expressions for the stress tensor and the charged current of 
the superfluid, as a function of the eight parameters characterizing flow. These eight parameters 
have three associated scalars, namely T, /i, ^. Symmetry considerations immediately allow us to write 
expressions for the stress tensor and current in terms of a set of arbitrary functions of these scalar 
fields. 

T^'^ = (p„ + p)u^u'' + p^r + i^ee 

= (p„ + P)u^'u'' + Pr;'^"^ + 
r = q,,u^^ - q,^j (2-1) 



The third equation in (2.1) is simply a statement of definitions. 

The reader may wonder why we have not allowed a u^^"^ cross term in the first of ( [2.l[ ). The 
reason is as follows. Were such a cross term to infact appear in the expansion, we could get rid of it 
by a redefinition of u'^. The assertion that no such cross term appears in the expansion of the stress 
tensor constitutes our definition of the normal velocity u^. 

The quantities p„, P, ps, Qn, Qs are all as yet arbitrary functions of the three scalar quantities 
(T^fi,^). The claim of the Landau Tiza two fiuid theory is that all these quantities may be derived 
from a single function, the pressure of these solutions, P = P{T, /i, ^) via the formulas 

Pn + P = QnP + Ts 
Ps = 

Ms - e - (2-2) 

dP = sdT + qsdfis + qndp, 
= sdT + qsdS, + qndp, 



-7- 



2.2 Relativistic Superfluid Dynamics 

At long distances (compared to an effective mean free path) a superfluid admits a fluid dynamical 
description. In this subsection we will present what appears to us to be the simplest and most logically 
satisfactory formulation of the theory of dissipative super fluid dynamics. Our presentation differs by 
a field redefinition from the more 'traditional' presentations of, for instance, Clark, Putterman and 
Landau Lifshitz 26, We will later describe the field redefinitions that allow us to transform 



our formulation to the more traditional one. 

In the relativistic context we choose the variables of superfluid dynamics simply to be u'^{x), the 
four velocity of the 'normal' fluid, T{x) the local temperature of the fluid, and £,'^{x), the local value 
of the gradient of the superfluid phase. More speciflcally, ^'^ is given in terms of the superfluid phase 
ip{x) by 

^/^(x) = (2.3) 

We have a total of eight fluid dynamical flelds. We will also sometimes use the terminology 

^l{x) = u{x)4{x). (2.4) 

We emphasize, though, that within the formulation presented in this subsection, fJ.{x) is not an in- 
dependent dynamical field, but merely terminology for the projection of ^(x) on the normal velocity 
field u{x). 

The equations of super fluid dynamics consist of the following relations 

a^T^" = 

J'' = (2.5) 
df^^i, - d^^i_, = 

The flrst two of these equations are simply the statement of the conservation of the stress tensor and 
the conservation of the charge current; these equations are true in any fleld theory. The last of these 
equations asserts that is the gradient of a scalar. In order to see that we have as many independent 



equations as variables, we note that (E^) is a local solution of the last equation in (2.5). This leaves 



us with 5 variables u^{x), T{x) and subject to the 5 remaining equations in the flrst two lines of 



While f^t(x) may be traded for ip{x) for counting purposes, the equations of fluid dynamics are 
more conveniently formulated directly in terms of the variables £,fj.{x) rather than the phase fleld ip^x). 
The reason for this is that gradients of "tpix) are not necessarily small in the regime of validity of 
superfluid dynamics, while gradients of ^'^{x) necessarily are. The introduction of the variables (,^{x) 
allows us to formulate the equations of superfluid dynamics in a systematic derivative expansion of all 
its participating flelds. 

The equations of superfluid dynamics constitute a closed dynamical system once they are combined 
with constitutive relations that determine the stress tensor and charge current in terms of the variables 
of fluid dynamics u'^,r, The constitutive relations take the form 

Ps 

e 



T^" = {Pn + + Plf" + 



(2.6) 



where tt'^'^ and J^^^g are respectively tensors and vectors that are first or higher order in an expansion 
in derivatives (of the fluid dynamical flelds) and all other quantities were deflned in the previous 
subsection. 



- 8 - 



2.2.1 A canonical fluid frame 

The equations of superfluid dynamics change their detailed form under field redefinitions. The quantity 
(x) has a microscopic definition. In the rest of this paper (x) will always refer to this microscopi- 
cally unambiguous value, and we will not permit arbitrary field redefinitions of ^^(x). 

On the other hand the fluid variables T{x) and u'^{x) are only really defined in thermodynamical 
equilibrium, and so allow for possible field redefinitions at derivative order (as derivatives parameterize 
departures from thermodynamical equilibrium). 

In order to completely specify the equations of superfluid dynamics we need to specify unambiguous 
definitions of the thermodynamical fields T[x), and u^{x). We specify these definitions by prescribing 
certain conditions on the dissipative terms tt^'' and J^^^^. We must choose these conditions so that 
they are not automatic, but can always be reached by an appropriate field redefinition, and completely 
fix field redefinition ambiguity. For instance we could work with the superfluid analogue of the 'Landau 
Frame' of ordinary fluid dynamics 

TT^."^ - (2.7) 

These 4 conditions serve to provide unambiguous definitions for the velocity, chemical potential * and 
temperature fields. For reasons that will become clearer below, we will refer to this choice of fluid field 
variables as the ^diss — frame. 

The equation (^^) may equivalently be written as 

T'^^M, = -Pnu^' + {u.Us)psK (2-8) 

where it is important that the functions p„ and ps on the RHS of (E^) are not independent functions, 
but are related to each other by the thermodynamical relations (|2.2|) ^ . 

2.3 'Fluid Frames' 

While the frame presented in the previous subsection seems to us to be rather natural from several 
points of view, it turns out not to be the fluid frame most commonly employed in analysis of superfluid 
dynamics. In this subsection we will describe a generalized framework for superfluid dynamics. In the 
next subsection we will describe how to transform between fluid frames. 

Conventional descriptions of superfluid dynamics are presented in terms of 9 fields subject to a 
single constraint rather than 8 independent fields. In this subsection we will take these 9 flelds to be 
T{x), p{x), M^(x) and ^''(a;)- 

While in the the previous subsection the fleld p{x) was simply convenient notation for u^(x)^^(a;), 
in this subsection fJ.{x) is an independent fleld variable; the relation between p{x) and u{x).^{x) is 
taken to be given by the so called Josephson equation 

u{x).^{x) = p{x) + Pdissix) (2.9) 

where i^diss (x) is a function of derivatives of fluid variables (i.e. it vanishes when all fluid variables are 
constants in spacetime). The quantity pdiss will be chosen in order to ensure that another condition 



*Note that here lJ,{x) is not an independent field and is given by(2.4). 

^In our presentation above we have used the temperature field as one of the independent fields of the superfluid 
dynamical description. Of course this choice is arbitrary; we could as well use any other thermodynamical field (for 
instance the energy density pn instead of the temperature, where pn is defined as the thermodynamical function of T, 
u.^ and 5) in place of the temperature. As such thermodynamical reparameterizations are ultralocal (i.e. do not involve 
derivatives) they do not affect the split of the stress tensor and current into perfect fluid and dissipative parts, and so 
do not constitute a change of frame. 
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we specify below is satisfied. Comparing ( 2A ) with ( |2.9| ) explains why we referred to the frame of the 
previous section as the ^diss = frame. 

The stress tensor and charge current continue to be given by the form (|]^), with the understand- 
ing, however, that all thermodynmical functions in those formulas are to be regarded as functions of 
the fields {T{x), iJ,{x),^{x)) rather than the fields {T{x),u.^{x),^{x)) (these two choices were identical 
in the previous subsection). As is not equal to u.^{x) for a general fluid flow, it follows that the 
perfect fluid part of the fluid stress tensor (at a given spacetime point) is, in general, not equal to 
stress tensor for any equilibrium flow (this is a complication that was absent in our formulation of the 
previous subsection). 

In the current formulation we have 5 thermodynamical fields - T{x), u'^{x) and fJ.{x) - that require 
precise definitions. This requires us to specify five equations (generalizing the four conditions in, e.g. 



(2J)) to give meaning to these fields. One natural choice ||29[ for these equations is 

We refer to the fluid frame defined by these equations as the Transverse Frame. These equations 
effectively determine the previously undetermined quantity Hdiss- 



Of course (2.10) is only one set of the possible set of five equations we must impose on our 



thermodynamical variables. Many other choices are possible. One other possible choice is 

f^diss — 0. 

which defines the fidiss = frame of the previous subsection. It follows that the formulation of 
superfluid dynamics, described in this subsection, is a generalization of the formulation of the previous 
subsection, and includes the later as a special case. 

We end this subsection by emphasizing our terminology. We refer to the formulation of fluid 
dynamics, presented in this subsection, as the formulation in terms of fluid frames. The key feature 
of this description of fluid dynamics is that the full microscopically defined gradient of phase, f , is 
taken as one of the variables of description. The perfect fluid part of the stress tensor and charge 
current is written in terms of and thermodynamical functions of f . Superfluid dynamics in a fluid 



frame is to be contrasted with super fluid dynamics in a modifled phase frames, introduced in §2.6 



2.4 Transforming between fluid frames 

In this subsection we supply the equations that allow us to transform between fluid frames. Let us 
suppose we want to make a change of variables that will take us from a completely unspecified fluid 
frame labeled 'our' to another frame labeled 'there', where the frame 'there' is a well deflncd fluid 
frame (e.g. the transverse or the iidiss = frame). We set 

'^there ~ '^our + 

Tthere = ^otir + (2.12) 
l^there — l^our ~^ 

The quantities 5u^, 5T and 5^ are necessarily of first or higher order in derivatives. As we work 
only to first order in derivatives in this section, we will effectively work only to linear order in these 
variations. 
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and M^^iir both be unit normalized, it is necessary that 
Suf.u'' = 

(at the order at which we work 5u^u'^^^ = 5u^u^^^^^ as the two differ at quadratic order in 5. Whenever 
an equation is true both of - say - uf^^^ and wj^/jg^e ^'^^^ simply omit the subscript in this subsection) . 
With infinitesimal variations restricted as above we find 



where by 



and 



OTT — n^j^j. T^there 



^f-^diss — {f^diss^our {f^diss^there 

and in these equations the symbol df(fi,^,T) represents the change in the function / under the first 
order variable change ( ^.12| ). These equations may then be used to obtain the 5 equations that 
determine the five unknowns J/z, ST and 6u^. For instance, if we are interested in transforming to the 
transverse frame (i.e. we want to set 'there' to be the transverse frame) we would require that T^tl^^^.^ 
and Jj'^g^g be orthogonal to u^f^^^^, giving the equations 

dpn + -^dps - vr^^^u^u,. 



(P + pn)5u^^ = ^dp^e - W'dpn - nZrU, (2.14) 



dqn + dqs- + {Jdiss)ourUf, = 



The 5 equations ( ^.14 ) determine Su^, Sp, ST, and so the change in the dissipative part of the 
stress tensor and current etc from (2.13). A similar procedure may be employed to transform into the 
Pdiss — frame or any other frame of interest. 

2.5 A canonical Entropy Current 

In this section we will define a 'canonical' entropy current in any fluid frame. Our definition is 

J's - su'^ ~ ^ JL. - ^ (2.15) 

where s is the thermodynamical entropy density of our fluid. 

Although this is not obvious, we have shown in Appendix ^ that this current is frame invariant. 
This means, for instance, that (2.15) deflnes the same vector fleld in the transverse as well as the 
Pdiss = frames. 

In the same Appendix we have also demonstrated that the divergence of this current is given, in 
any fluid frame, by 







7i" 


.Y. 




.T. 



JLs + ) (2.16) 
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2.6 Modified Phase frames 



As we have emphasized above, the theory of superfluid dynamics formulated in any fluid frame in 
which ^diss 7^ (like the transverse frame) has slightly hybrid features. The perfect fluid part of the 
charge current and stress tensor involves thermodynamical functions of /i not u.^, but is written (in 
index structure) in terms of the full field This ensures that, for a general fluid flow, the prefect 
fluid stress tensor and charge current, at any given point, does not equal the stress tensor and charge 
current for any equilibrium flow. 

Fluid dynamics formulated in a modifled phase frame eliminates this slightly unpleasant feature 
by working with a modified gradient of phase field, ^q{x) defined as 

= -fiuf" + wf" 

^ (2.17) 

Co = V /^^ - ^2 

where is defined as the part of projected orthogonal to . 

Note that is not equal to the phase field (because /i 7^ u^^^ ). Instead the relationship 
between these two fields is given by 

, , ^ M (2-18) 

Note also that, by construction, ^q.u = fi. 

The fiuid stress tensor, charge current, in a modified phase frame, are assumed to take the form 



T^" = {Pn + PW'W' + Pr;"'^ + gCo^Co + 

so 

t.0 



(2.19) 



where all thermodynamical functions are taken to be functions of (T, /i,Co)- 

As in the case of fluid frames, the precise definition of any given modified phase frame requires the 
specification of 5 additional conditions (to give precise meaning to the fields T{x), u^{x) and /i(x)). The 
'Landau-Lifshitz-Clark-Putterman' frame is a modified phase frame in which the additional conditions 
are taken to be 

As the name suggests, this is the frame employed by Landau Lifshitz, Clark and Putterman ||2^, ^ 
^(in a non relativistic context) in their analysis of superfluid dynamics. 



It is not difficult to generalize the analysis of §2.4 to describe the transformation from a modified 



phase frame to a fluid frame or vice versa. We will implement such a frame change in the next 
subsection. 



2.7 The canonical entropy current in modified phase frames 

The most general modified phase frame may be obtained starting from the most general fluid frame 
and then reexpressing all quantities in terms of Co a-nd Co rather than C^ and C- Let the fluid frame we 
start with be characterized by tt^"^ and J^iss ^'^'i Pdiss- The generalized phase frame, obtained from 
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this fluid frame by reexpressing i^'^ as a function of has 

JLs = JLs - dqnu^ + die, + I fid^ssu^ ^2 20) 



ISS 



where / — = Further, any thermodynamical function A in the fluid frame is related to the 
corresponding thermodynamical function A in the modified phase frame by 

dA = A-A = Ai^o) - AiO = -^^d^ss (f^ 



As we have explained above, the canonical entropy current in an arbitrary fluid frame is given by 
( ^.15 ). Applying the transformation formulas described above, this entropy current may be expressed 



in terms of the modified phase thermodynamical and dissipative parameters as (see Appendix A.3.1) 

'J-^ " diss rji • rp 



In obtaining ( 2.21 ) we have used the thermodynamical identity ( p.2| 



It is also possible to transform the equation for the divergence of the entropy current, ( 2.16 ), into 



the modified phase frame (or simply rederive this expression directly in the modified phase frame). 



We find (see Appendix A. 3) 







7^" 


.Y. 




.T. 



JLs+f^d^ssP^''d,(l^] (2.22) 



3. A theory of first order dissipative superfluid dynamics in fluid frames 

In this subsection we will present a 'theory' of dissipative superfluid dynamics to first order in the 
derivative expansion. By this we mean that we will parameterize the allowed forms of tt'^'^ and Jj^.^^^ 
at first order in the derivative expansion. Our parameterization will be in terms of a certain number 
of undetermined functions of T /i and ^. One of these functions is the viscosity of the normal part of 
the superfluid. Following standard (but slightly misleading) usage, we will refer to these functions as 
dissipative parameters of the superfluid. 

3.1 Summary of arguments and results 

As the analysis of this section will be rather lengthy, we flrst present a summary of our logic and our 
procedure. To start with we simply classify all onshell inequivalent one derivative contributions to 
TT^'', J^jjjj, and fj^diss- It is not difficult to establish that, in any given frame (e.g. the transverse frame 
or fidiss = frame or the Landau-Lifshitz-Clark-Putterman frame) there exists a 36 parameter space 
of inequivalent first derivative corrections to the equations of superfluid dynamics (assuming parity 
invariance) . 

In order to cut down the set of possibilities we then follow Landau Lifshitz, Clark and Putterman 



||25[ g6|, |27| to make the central assumption of this subsection. We assume that the entropy current 
takes the canonical form described in §2.5 ^. A local form of the second law of thermodynamics then 



^This assumption is not universally valid. In particular it seems certain to fail in situations in which the U{1) 
symmetry in question has a U{1)^ anomaly. It may also fail in other circumstances. We leave the investigation of the 
validity of this assumption to future work. 
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asserts that the equations of superfluid dynamics should be geared to ensure that the divergence of this 



entropy current is positive. Using the expression (2.16) we find that this requirement cuts down the 36 



parameter space of possible one derivative corrections to the entropy current to a 21 parameter space 
of possibilities. The coefficients in this 21 parameter space are further constrained by a complicated 
set of inequalities that ensure positivity of the entropy production. (One of these inequalities, for 
instance, asserts the positivity of the normal viscosity). Finally, the Onsager reciprocity relations 
relate 7 of the remaining parameters to each other, leaving us with a 14 parameter space of dissipative 
coefficients. As mentioned above, these 14 parameters (each of which is a function of T, fj, and ^) 
are further constrained to obey a set of inequalities. As far as we are aware, there are no further 
restrictions on this 14 parameter space from general principles. 

To end this summary we explain how the framework presented in this subsection relates to previous 
work. The programme outlined in the paragraph above was implemented by Landau and Lifshitz [ p7[ 
for the special case of flows with zero (or negligibly small) superfluid velocities. Landau and Lifshitz 
found a set of equations with 5 first order dissipative parameters. Our 14 parameter set of equations 
indeed reduce to the Landau Lifshitz form upon setting the superfluid velocity to equal the normal 
velocity; consequently our framework agrees with that of Landau and Lifshitz within its domain of 
validity. 

Clark and Putterman ^ extended the Landau Lifshitz programme to the case of nonzero 
superfluid velocities. The end result of their analysis was a thirteen parameter set of equations. We 
believe that Clark and Putterman overlooked one allowed parameter 

Reinstating that parameter yields our 14 parameter set of equations. Thus Clark and Putterman's 
equations are a special case of our 14 parameter equations with one parameter set to zero. 

As we have explained in the previous section, one of the complicating features of superfluid dy- 
namics is that one can work in many different frames. In this section we work out the 'theory' of 
dissipative superfluid dynamics in the two natural fluid frames described in the previous subsection, 
namely the transverse frame and the ^diss — frame. In the next section we present the equivalent 
analysis for in the Landau-Lifshitz-Clark-Putterman frame. 

3.2 Counting of parameters 

At a given spacetime point, a superfluid flow has two independent velocities; the normal velocity 
and the superfluid velocity These two velocities break the local Lorentz rotation group SO {3, 1) 
down to 5*0(2), the group of spatial rotations in the plane orthogonal to both velocities. 

In § |7.2| we enumerate the onshell inequivalent first derivatives of all fluid dynamical flelds at a 
point. We found it convenient to classify these derivatives as scalars (spin 0), vectors (spin ±1) and 
tensors (spin ±2) under the unbroken 5*0(2) described in the previous paragraph. As explained in 



§7.2 it turns out that there are 6 onshell inequivalent parity even scalars, 5 onshell inequivalent parity 
even vector and 2 onshell inequivalent parity even tensor first derivatives of fluid dynamical fields. ® 
In order to be specific we will assume in the rest of this subsection that we are working in the 
transverse frame. Very similar arguments can be made in a fluid or modified phase frame, and give 
identical results. 



^Specifically, the traceless symmetric 3 index tensor listed in equation (A VI-9) of Putterman's book is not unique. 
Another such tensor is given by 

Yijk = Wi{wjWk - (l/3)io^(5jfc). 

*In addition we have one additional parity odd scalar field. Further, every vector can be transformed to a psuedo 
vector Vfi according to the formula Vfi = e^j^^^y n"u''. 
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In the transverse frame, tt^'' has two inequivalent scalar components C^^'^Tr^jy and tt^, one vector 
component P^Tr"'^^^ and a single tensor component Pai.iPpu'^"^ , where Pq/? is the projector orthogonal 
to the u^, plane. On the other hand J^j^^ has one scalar component ^^j^^C/i and one vector 
component Pafj.Jdiss- Finally /idiss has one scalar component. It follows that the total number of 
undetermined parameters in the arbitrary expansion of tt^", J^j^g and f^diss in terms of first derivatives 
of the fluid dynamical fields (assuming parity invariance) is given by 

4x6 + 2x5 + 2 = 36 

where the three terms above originate in the scalar, vector and tensor sector respectively. ^ 

3.3 Constraints from positivity of entropy production and Onsager relations in the trans- 
verse frame 

In this subsection we will explore the constraints on dissipative coefficients from the physical require- 
ments of positivity of entropy production and the Onsager reciprocity relations. We will find these 
requirements cut down the 36 parameter set of possible dissipative coefficients (assuming parity in- 
variance) to a 14 parameter set of coefficients that are further constrained by positivity requirements. 

For concreteness we present our analysis in the transverse frame. We will record the results of 
similar analysis in other fluid frames in later subsections. 

3.3.1 Constraints from positivity of entropy production 



The divergence of the 'canonical' entropy current, given by ( 2.16 ), involves only terms proportional to 
d^u^i:^'^ , d,^{^/T)J^^^^ and ^idissd^{qs£,^ / C)- Let us examine these terms one by one. In the transverse 
frame 



where cr^^ is the traceless symmetric part of d^u^, projected in the direction perpendicular to u'^ . 

_ papP f dgup + dpug [d.u] 

and ' (3-1) 

P^^u = The projector = rj^^, + u^u^, 

Now the field fi^ti/ has one scalar piece of data 



one vector piece of data 
and a tensor piece of data 



' p.f u u ^ a 







^Dropping the assumpti on of parity invariance we have 4x7 + 2x10 + 2 = 50 independent dissipative coefficients, 
^"in the terminology of §7 

y below, [Vt]^ = ^[Vs], 



^^In the terminology of S 
^•^In the terminology of § 



2 L 

7.2 below, = ^[Ti]^^ 
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The trace of n'^" couples to another scalar piece of data 



13 



In the expressions above n/^ is the unique normal vector in the plane spanned by and that 
is orthogonal to and is given by 



w 



with w being the norm of the vector. 
Similarly, in the transverse gauge 



where P^v is the projection operator (defined in (3.1)) that projects orthogonal to it^ only . The 
quantity P^d^{ii/T) has one scalar piece of data 



and one vector piece of data 



Sb - {n^d^) ifi/T) 



Finally 



Sa 



T3 



is itself a scalar piece of data. 



In other words we conclude that the expression for the divergence of the entropy current, (2.16), 
depends exphcitly (i.e. apart from the dependence of tt'^" J^-^^ and Hdiss on these terms) only on 4 
scalar expressions, 2 vector expressions and one tensor expression. Let us choose these 4 vectors scalars 
Sa, Sb, Sw and S^', supplemented by 3 other arbitrarily chosen scalar expressions S'^ (m = 1 ... 3) as 
our 7 independent scalar expressions. Similarly we choose the 2 vectors [Va]^ and [Vf,]^ supplemented 
by 3 other arbitrarily chosen expressions {m = 1 ... 3) as our four independent vector expressions. 

We also choose Tf^i, as one of our two independent tensor expressions We proceed to express tt^''', 
"^diss ^'^'^ Mdiss as the most general linear combinations of all combinations of independent expressions 
allowed by symmetry 



^^In the terminology of §7.2 below, S^i = Sa + Si^. 

^"^We could now go ahead and use the perfect fluid equations to solve for for the dependent data in terms of independent 



data; however we will not need the explicit form of this solution in this subsection 



- 16 - 



rp6 



PaSa + PbSb + P-ujSw + Pw'Sw' + ^ ' Pm^rr 



P 



rn—l 
3 

' / ^ m m 
m — 1 



'2 



(3.2) 



diss 




RbSb 



Pw Su; 
3 



CbVi: 



m— 1 



QaSa + QbSb + QwSw + Qw'Sw' + ^ ^ Qm^ra 



Plugging (|3l8|) into ( p^ we now obtain an explicit expression for the divergence of the entropy 
current as a quadratic form in first derivative independent data. We wish to enforce the condition 
that this quadratic form is positive definite. Now the quadratic form from ( 2.16| ) clearly has no terms 
proportional to (5*4,)^. It does, however, have terms of the form (for instance) SaS^, and also terms 
proportional to S^. Now it follows from a moments consideration that no quadratic form of this 
general structure can be positive unless the coefficient of the SaS^^ term vanishes. Using similar 
reasoning we can immediately conclude that the positive definiteness of ( 2.1(\ ) requires that 



EL 



Cm — Rm — T2 — 0. 



(3.3) 



(3.3) is the most important conclusion of this subsubsection. It tells us that a 21 parameter set of first 
derivative corrections to the constitutive relations are consistent with the positivity of the canonical 
entropy current. 

Of course the remaining 21 parameters are not themselves arbitrary, but are constrained to obey 
inequalities in order to ensure positivity. In order to derive these conditions we plug (3.3) into (3.2) and 
use (2.16) so that the divergence of the entropy current is the linear sum of three different quadratic 
forms (involving the tensor terms, vector terms and scalar terms respectively) 



d^J^ = (g, + Qv+ Qt) 



(3.4) 



where 



Ca 



Va 



Qt - -tT^ 

5a) T 

Cb + Ea 



CaVl - {Cb + Ea) VbVa ~ E^V^ 



2Ca 



Vb 



Eb 



{Cb + Eg 

^Ca 



(3.5) 



^^For instance the quadratic form + cxy (where c is a constant) can be made negative by taking ^ to either positive 
or negative infinity (depending on the sign of c) unless c = 0. 
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Qs 



PwS^ Tr^^'S^i Qa^a P-b^b 



— {Qw + Pa) SwSa ~ iQw' + To) SyjiSa — {Rw + Pb) S^Sb (3.6) 

— {Rw' + T{) Sw'Sb — {Ra + Qb) SaSb — {Tw + Pw') SwSw' 

Positivity of the entropy current clearly requires that Qt Qv and Qs are separately positive. Let us 
examine these conditions one at a time. For Qt to be positive it is necessary and sufficient that r < 0. 
This is simply the requirement that the normal component of our superfluid have a positive viscosity. 
In order that Qv be positive, it is necessary and sufficient that 

Ca<0, Eb<0 and > {Cb + Eaf . (3.7) 

Note that this expression involves Ca and Eb on the LHS but the different quantities Cb and Ea on 
the RHS; the last inequality above is satisfied roughly, when Cb and Ea are larger in modulus than 
Ca and Eb- 

Finally Qs, listed in (B.6), is a quadratic form in the the 4 variables Sa, Sb, Sw and Sw'- We 
demand that this scalar form be positive. We will not pause here to explicate the precise inequalities 
that this condition imposes on the coefficients. See below, however, for the special case of a Weyl 
invariant fluid. 

3.3.2 Constraints from the Onsager Relations 

In the previous subsection we found that flrst order dissipative corrections to the equations of perfect 
superfluid dynamics take the form 



TP 

diss 



{PaSa + PbSb + PwSiu + Pw' Sw') { n^Ui, — 
+ {TaSa + TbSb + T.fjuSw + Tw'Sw') P'^" 

+ Ea {Va^n" + V:n^) + E^ {V.^'n'' + Vj^ n^) 

{RaSa + RbSb + RwSw + Rw'Sw') 



(3.8) 



+ CaV^ + CbV^ 



fJ-di. 



' Su 



[QaSa + QbSb + QwSw + Qi> 

where the coefflcients in these equations are constrained by the inequalities listed in the previous 
subsubsection. The coefficients that appear in these equations are further constrained by the Onsager 
reciprocity relations (see, for instance, the text book for a discussion). These relations assert, 
in the present context, that we should equate any two dissipative parameters that multiply the same 
terms in the formulas ( p.5| ) and (3.6) for entropy production. This implies that 

Qia — Pa 7 Qio' — Pai Rw Pb 
Rw' Pb, Ra — Qb} Pw — Pw 

and 

Cb — Ea 

In summary we are left with a 14 parameter set of equations of first order dissipative superfluid 
dynamics. The requirement of positivity constrains further these coefficients to obey the inequalities 
spelt out in the previous subsubsection. 



(3.9) 
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3.4 Weyl Invariant Superfluid Dynamics in the transverse frame 



Let us now specialize the results of §3.3.1 and §3.3.2 to the case of super fluid dynamics for a conformal 



superfluid. The analysis of §3.3.1 is simplified in this special case by the fact that the trace of the stress 
tensor vanishes in an arbitrary state (and so in the fluid limit) of a conformal field theory. This fact 



reduces the number of explicit scalars that appear in (2.16) from 4 to 3 (the scalar Sw' never makes an 
appearance). It follows that the requirement of Weyl invariance forces P^j' = Rw' — Tw — Q^,i — 0. 
Moreover the requirement that tt^'' be traceless forces Ta = T\j = = 0. It turns out that there are 
no further constraints from the requirement of Weyl invariance. The expansion of the dissipative part 
of the stress tensor and charge current for a conformal superfluid is given by 



[PaSa + PbSb + PwSw) (^nf^Ui, - 
(RaSa + RbSb + RwSw) 



+ CaV^ + CbV^ 
[i'disa = ~ [Qa'S'a + QbSb + QwSw] 

The entropy production is given by 



where 



Qt = -tT^ 

QV=- CaV^ - {Cb + Ea) VbVa ~ EbV^ 



= -Ca 



Va 



I Cb + Ea 



Vb 



Eb 



{Cb+Eaf 
^Ca 



(3.10) 



(3.11) 



V 2a 

Qs = ^ PwS^ — QaSa — RbS^ 

+ (Qui + Pa) SwSa — {Ra + Qb) SaSb + {Rw + A) S^Sb 

For the entropy current to be positive it is necessary and sufficient that t < and that 
a < 0, Eb<Q and mCa > (Cb + Ea)^ . 

and that the quadratic form 

Qs = aixl + 022^2 + + biXiX2 + b2X2X3 + baxixs 

-\ 2 



(3.12) 
(3.13) 

(3.14) 



ai 



bi 



as 



4a 1 



2^3 



2a, ' 

2ai62 - bibs 
4aia3 — 5? 



(3.15) 



(4aia2 - 6?)(4aia3 - bj) - (2ai62 - bib^)' 



4ai(4aia3 — 63) 
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is positive with xi = S^, X2 = Sa and = Sb and 

ai = -Puj, a2 = -Qa, as^-Rb, h = Qu, + Pa, b2 = -{Qb + Ra), h = + Pb 
For the last quadratic form to be positive it is necessary and sufficient that 

ai > 

4aia2 > hi (3.16) 
(4aia2 - fo?)(4aia3 - hi) > {2aib2 - ^163)^ 

By rewriting ( ^.151 ) as a sum of squares in a cyclically permuted manner we can also derive the cyclical 
permutations of these equations. 

In summary, the most general Weyl invariant fluid dynamics consistent with positivity on the 
entropy current is parameterized by a negative ri, 4 parameters in the vector sector constrained by 
the inequalities ( |3.14 ) and 9 parameters in the scalar sector, subject to the inequalities ( |3.16 ). These 



(3.17) 



14 dissipative parameters are further constrained by the 4 Onsager relations 

Qw — Pa, Rw — Pb, Ra — Qb, Cb — Ea 

leaving us with a 10 parameter set of final equations. 

3.5 Weyl Invariant Super Fluid dynamics in the ^diss = frame 

In this section we present constitutive relations of super fluid dynamics in the frame that we consider 
the most natural, namely the ^diss — fluid frame. We only present our final results, and specialize 
to the case of a Weyl invariant fluid for simplicity. We find 



[PsS + PbSb + Pn,S^] 



TP- 

diss 



+ Ea {Vl^Tf + V:n^) + Eb iy^n- + V^n^^) 

[QsS + QbSb + QwS^] 

+ [RsS + RbSb + Rn^S^,]nt' 



(3.18) 



+ CaV,^ + CbV^ 

where S = (u.d) (^) 

The equation of entropy production is given by 



d^J^ = (Q, + Qv + Qt) 



(3.19) 



where 



}V^- CaV^ - {Cb + Ea) VbVa - EbV^ 



Ca 



Va 



Qt 

Cb + Ea 
2Ca 



Vb 



Eb 



jCb+Ea)' 
iCa 



Vb' 



(3.20) 
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and 



Qs — — PwSl^, — QgS^ — RbSl 

+ (Q^ + Ps) S^S + (i?^ + Pb) S^Sb - {Rs + Qb) SSb 



(3.21) 



The positivity of Qt is equivalent to the requirement that ri < 0. The positivity of Qy is equivalent 
to the condition 

Ca<0, Eb<0 and AEbCa > (Cb + E^f 



The implications of the requirement of the positivity of Qs are exactly as in §3.4. The Onsagar 
relations imply that 

Qw = Ps, Rw = Pb, Rs = Qb, Cb = Ea 

Once again we have a 10 parameter set of dissipative equations of super fluid dynamics. 

4. Dissipative Superfluid dynamics in the Landau-Lifshitz-Clark-Putterman 
Frame 



In the previous section we have presented a 'theory' of first order dissipative super fluid dynamics 
in the transverse fluid frame. In this subsection we present the equivalent 'theory' in the Landau- 
Lifshitz-Clark-Putterman modified phase frame. As the logic of this presentation mirrors that of the 
previous subsection in every way we will omit all derivations and present only our final results. After 
we have presented our results we will explicitly take the non relativistic limit and compare with the 
analogous expressions in Clark, Putterman and Landau and Lifshitz [|5, 26, 



4.1 Allovired forms for dissipative terms 

As we have explained above, the Landau-Lifshitz-Clark-Putterman frame is a modified phase frame 
that specified by the relations 



0, JLs = o 



We take the stress tensor to be given by 



{Qf^u" 



^"u'') + HP'"' 



where 



(4.1) 



(4.2) 



The divergence of the entropy current in this frame is easily obtained from (A. 20) and takes the 



form 



8 



[{d.u)n + nfVp.] Q" [duT + T{u.d)u,] 



T 



r2 



f 

T 



(4.3) 



The first derivative expressions of fiuid dynamical fields that appear explicitly in (O) may be 
written in terms of the following scalars 



Si = 71 — , Y.2 = n''af_,yn , = d^w^, 1:4 = 



d. 



(^) 



T 



T2 



the following vectors 



T 
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and the following tensors 
where 



w 

In terms of these expressions, positivity of the entropy current requires 



n = [ttiSi + n2^2 + TTaSs + ^4^4] 



(4.4) 



nr = (Pi Si + P2S2 + + P4S4) (n^n^ - ip"'^ 
+ El (Wf + VVi'n'*) + E2 {W^n" + W^n'') 



With this form for the dissipative parameters, the divergence of the entropy current is given by 

df,J^ = {Qs + Qv + Qt) 



where 



Qs = - gi'^S? - P2S2 - ^3^1 - I14T.I 

- {Q['^ + /Xi)SiS4 - (P4 + M2)S2S4 - (774 + ^3)^3^4 

- iQi''^ + Pl)SiS2 - {Qi'^ + 7ri)SiS3 - {W2 + P3)S3S2 



Qv = - Q^^'^Wi^ - (q^^) + El) W1W2 - E2W, 



(4.5) 



Wi 



W2 



E2 



4Qi 



(") 



The positivity of entropy current requires that the each of the three quadratic forms Qs, Qv and Qt 
are positive definite. 

The Onsager relations take the form 



Qi'^=Pl, Qi"^=7Tl, P3=W2 



Q 



2 

(v) 



(4.6) 



El 



Once again we have 14 indpendent dissipative coefficients, constrained by inequalities that follow from 
the requirement of positivity of the entropy current. 
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4.2 The Limit of zero superfluid velocity 

In this subsubsection we will examine how the description of the previous subsection must simplify in 
the limit — > 0. In this limit the superfluid and normal velocity are coincident. In this limit the vector 
n'^ no longer has any significance. The equations for the dissipative parts of the stress tensor, current 
and fidiss must all possess the full 5*0(3) invariance of rotations that leave m'" fixed. The expressions 
for all dissipative quantities must be analytic functions of the projected superfluid velocity vector field 
— Pj^S,^ — TcCn^- This requirement imposes several constraints on the C ^ behaviour of the 
coefficients in dissipative terms, including 

lim q["^ = lim g^'^ = hm qI"^ = lim Q^"' = 

w^O w~^0 10-^0 uf— ^0 

lim Pi = lim P3 = lim P4 = lim E2 — 

lim fii = lim /i2 — lim vri — lim 112 — (4 7) 

to— >0 lu— )-0 to— >0 to— >0 ^ ' 

lim P2 = lim Ei = lim t = rj 

tu— >0 lu— >0 lu— >0 

lim g^"^ = lim Qi"^ = K 

to— >0 to— J-O 

16 

It follows that, in this limit the 21 independent dissipative coefficients of the previous subsubsection 
(before imposing Onsager relations) reduce to only 6 non-zero coefficients 

?7, K, /^3, /^4, TTs, 7r4. 



The relations in (4^) simplify to 

^J'diss — + 



n = [^3^3 + 7r4E4] 



(4.8) 



The Onsagar relation set fi^ = tt^, finally leaving us with 5 dissipative coefficients, as predicted by 
Landau and Lifshitz [p7| . 

4.3 The Non Relativistic Limit 

In order to make contact with the results of Clark and Putterman |2^(and with those of Landau 
and Lifshitz in the limit of zero superfiuid velocity) we will now explicitly take the non relativistic 
limit of the formulas presented in §4.1. 

'^^To see where these conditions come from, consider, for example, a constitutive relationship for a quantity vi which 
we know is a full SO{3) vector when ^ vanishes, but is simply the sum of an 50(2) vector and an SO{2) scalar at 
nonzero C^. At arbitrary values of f , the scalar and vector parts of vi can be expanded as arbitrary linear combinations 
of all independent SO(2) scalars and vectors. However as — > 0, the scalar component can be expanded in only those 
scalars that arise from the decomposition of an 50(3) vector. As all 50(2) scalars are not of this form, this implies 
the vanishing of a number of coefficients in the expansion. Moreover the coefficient of these scalars has to equal the 
coefficient of the corresponding SO (2) vectors in the vector part of vi. This requires that two otherwise independent 
coefficients are equal. Similar remarks hold for the expansion of vector fields. 
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The non relativistic limit is taken as follows. We set 



{l,v}, w^ = {0,w}, n'^ = {0,n}= 0,- 



w 

(we ignore all terms of quadratic or highe order in v and lu) and 

Note in particular that all vectors and tensor projected orthogonal to in this limit, are purely 
spatial. 



The various scalar, vector and the tensor terms listed in §4.1 reduce to the following expressions 
in the non relativistic limit: 



Si = 



{n.d)T 
T ' 



f w 
T 



diT 

[Wi]-* = [^^1* " ^ "'^^ 



(4.9) 



(note that all vectors and tensors are purely spatial, and we have only listed their spatial components) . 
It follows that 



(4.10) 



where all the quantities in the equation above are given by plugging (4^) into (4^). The resultant 
equations are very close to equation AVTll - AVT14 in Putterman's book The main difference 
between our equations and those of Putterman is that our equations have 21 free parameters, in 
contrast with Putterman's 20 free parameters. The equations of this subsection reduce to (A VI -11 



to A VI-14) in [g6| if we set Q^' 



(^) 



0. 



The Onsager relations of § [4.1| carry over unchanged to the non relativisitc limit. These relations 
are identical to the Onasager relations imposed by Putterman |2^ and reduce the number of dissipative 
coefficients to 14 in our formulation, but to 13 according to Putterman [ p6| . 

The non relativistic version of the tZ; — > constitutive relations (4.8) is simply given by by plugging 



( [4.9| ) into ( |4.8| ); the resulting 6 parameter set of constitutive (cut down to a 5 parameter set by the 
Onsager relations) agree with those presented by Landau and Lifshitz ||27|] . 



5. Equilibrium Superfluid Thermodynamics from Gravity 



In this section we explain that the equilibrium superfluid thermodynamics, reviewed in §2.1, may be 
derived from gravity very simply. This section, which is abstract in nature, is a generalization and 
reworking of the beautiful paper of Sonner and Withers M on the same subject. 
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Consider a general gravitational system governed by the action 

S = Y^ J V5(^?+12-^(V^iW*)Fa6i^'^'' + V2(0<^*)D,</.5'^0* + V3W*))) (5.1) 

where = V — iAa and I)q = V + iAa and all the potentials Vi, V2, V3 are real. We assume 
that this system admits a homogeneous stationary asymptotically AdS family of solutions - dual to 
homogeneous stationary superfluid flows - that take the form 

Metric : ds'^ = -2g{—)Ufj^dx'^dr - rlf{ — )u^j,u^dx^dx'' + k{—)n^nudx^dx'' 



Tc rc rr 

+ rl j{ — ) [rinUi, + UnTiy) dx^dx^ + Pn^dx^dx'^ , 
rc 



Gauge field : r^A = H( — ) u'^d^ + U — ) n^'^. 

rc rc 

r 
rc 



(5.2) 



Bulk scalar field = 
where 

Pnv = r]^v +u^Uv - n^riu. (5.3) 
Here and are two arbitrary constant vectors obeying 

u^n^ = 0; u^Uf, = -1; 71^71^ = 1. (5.4) 

We work in a gauge such that the scalar field is real 4>* — cj) (this implies A"^ — 0), so that the 
boundary value of the gauge field gives the superfluid velocity. We choose the constant vector so 
as to ensure that the killing vector coincides with the generators of the event horizon of our solution. 
TT,^ is then uniquely determined by (5.4) together with the requirement that A^ at infinity (i.e. ^^) 
can be written as a linear combination of u'^ and . 

We now choose coordinates so that the killing vector points along the direction of the and 
the vector dx points in the direction of n'^ . Our solution retains rotational invariance in the remaining 
two spatial directions. We also work in the rescaled variables f- and r^a;^ in terms of which ( ^^ ) 
reduces to 

ds^ = 2.9 (r) dv dr - f{r)dv'^ - 2j(r) dv dx + k(r) dx^ + dyl 

A'- = 0, A" = H{r), = L{r), A^' = 0, A' = Q (5-5) 
Bulk scalar field = (/)(r) 

The solutions above are translationally invariant in the field theory spacetime directions. They are 
characterized by an onshell action density (action (^^) per unit volume), a stress tensor and a charge 
current. Merely from symmetry, and with an appropriate definition of , the stress tensor and the 
charge current of our system necessarily take the form ( ^.1[ ). In Appendix ^ below we demonstrate 
that the functions that appear in (^^) obey all the thermodynamical relations of the Landau Tiza 
two fluid model described above. 

While the details of our demonstration are presented in Appendix ^ we very briefly review the 
main ideas here. 

Our system has three inequivalent translationally invariant killing vectors: u^d^ = d^, 
and fc''(9p where fc'^ is transverse to the normal as well as the superfluid velocity. As is standard in 
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derivations of black hole thermodynamics, we demonstrate using the equations of motion that 
(where 6'^ is any of the vectors above) takes a particularly simple form. The resultant identities allow 
us to demonstrate the onshell vanishing of three total derivatives built out of the functions in the 



metric and gauge field (5.2). Integrating these expressions allows us to prove the Smarr-Gibbs-Duhem 



relations listed in the first two of equations (2.2). We are also able to show that the killing vector 
that reduces to the null generator of the horizon of our solution also defines the normal velocity of our 
fluid. Recall the later is defined to ensure the absence of u^u'^ cross terms in the stress tensor. 

Next we follow [Q to rewrite the action of an equilibrium configuration of the sort studied in this 
section as the integral of a total derivative. Performing the integral we find that the action evaluates 
to the difference of two terms, one at infinity and the second at the horizon. The term at the horizon 
vanishes while the term at infinity evaluates to the negative of the pressure of the solution. 

Finally, the action density (hence the pressure) of this system may be regarded as a function of its 
temperature T, chemical potential /i and superfluid chemical potential fis. It is well known in classical 
mechanics that the onshell action as a function of initial and final coordinates, and obeys the 
relation dS = PodqQ — P^fdqy. Sonner and Withers have demonstrated that the application of this 
result to our system (with r being regarded as the time) gives 

dP = qndu + qsdiis + sdT (5.6) 

where s is the entropy density of the solution given by the Hawking- Beckenstein formula. 



6. An analytically tractable limit of hairy black branes 

As we have explained in the previous section, general gravitational methods are powerful enough to 
tightly constrain the gravitational solutions dual to fiuid dynamics, in particular to demonstrate that 
the thermodynamics of these solutions is given by the Landau-Tiza two fiuid model. In order to 
explicitly determine the thermodynamics of any particular gravitational system, however, we need 
to explicitly determine the solutions dual to uniform superfiuid fiows. Unfortunately, the ordinary 
differential equations that arise in this attempt have proved so complicated that it has not proved 
possible to analytically solve for hairy black branes (the gravitational duals to superfiuids) in any 
reasonable gravitational system The only analytic results that we are aware of, for hairy black 
brane solutions, are those of Herzog [^. Herzog considered a very special model, the model of a 
charged scalar field of = —4 and infinite charge e (i.e. a model in the so called probe limit). He 
demonstrated that this model displays a second order phase transition towards superfluidity whenever 
l#l>2. 

When 1^1 is just larger that 2, the stable gravitational solutions develop a scalar vev. Let e denote 
the value of this vev. Herzog was able to generate the relevant gravitational solutions perturbatively 
in e and also perturbatively in the difference between superfluid and normal velocities. 

In this paper we will be interested in probing the structure of viscous superfiuid dynamics from 
gravity. In the infinite charge or probe limit of scalar and gauge dynamics do not back react on 
spacetime. In order to probe the dynamics of the interaction between the stress tensor and the charge 
current we need to go beyond the infinite charge probe limit. In this section we generalize Herzog's 
perturbative construction of gravitational solutions to go beyond the probe limit. In other words we 
generalize Herzog's infinite e solutions to retain the first nontrivial correction in a (i) expansion. 

In the next section we will use the results of this section as an input into the fiuid gravity map, 
in order to generate gravitational solutions dual to viscous superfluid flows. 

^''Of course much attention has been focused on the numerical solutions of the relevant equations in several models. 
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6.1 The bulk system and the equations of motion 

Following Herzog [|| we consider the system 

^ = 16^ / ^'''^ + 12 + ^ I^-\f,'^F'"' - l\D,^\' + 2|0r)) , (6.1) 

Where 13^ — — iA^^, and is the gravitational covariant derivative. Note that this system is a 
special case of the system (5.1) (where we have chosen particular simple forms of the functions Vi,V2 
and V3). 

The equation of motion for the scalar field and the gauge field that follows from (3.1) are respec- 
tively 

D^D^c/) + 40 = 0, (6.2) 

and 

D^F'^" = Ij", (6.3) 
where the current = i {(j)*Dfi(f) — (t>{D ^(j))*) . The Einstein Equation that follows from ( |o| ) is 

G^u — 6.9^1/ = ^ ((Tmax)/ji/ + (7niat)/ji/) , (6-4) 



where 



(6.5) 



6.2 Boundary Conditions and Solutions 



We search for solutions of the form (5.5). The 4- vectors rt^, defined in the previous section, may be 
computed as follows. Let 

It follows that is given by 

As explained in the previous section we choose gauge so that our scalar field is real (j){r) = <jf{r) 
by a choice of gauge. This gauge choice, plus the static nature of our solution, forces A'" = (as can 
be seen by comparing the equation of motion of with that of <jf). 

We search for solutions that obey the following large r boundary conditions 



k{r) 




fc2 


fir) 


= rN 


_^ 




- ii 


f . . 


L{r) 




f . . 




e 


f . . 



(6.6) 
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It turns out that the conditions above, together with the equations of motion, automatically ensure 

lim g{r) = 1 

SO that this condition, while true, does not have to be additionally imposed. Also, it turns out that 
the coefficient of term in the asymptotic expansion of H{r) is fixed by equations of motion and 
the requirement that 4> be regular at the horizon. 

Our functions are also constrained at r = 1 as follows 



= /(I) = (6.7) 

On the other hand the functions H{r), k(r), L{r) and 0(r) are required only to be regular at r = 1. 

It is possible to argue that there exists an 8 parameter class of solutions of the form (pT^), to the 
system (pH]), subject to the boundary conditions listed above. One of these parameters is Tc in ( ^.2] ). 
The three normal velocity parameters can be set to zero by a boost, and rotations can be used to 
point the supcrfluid velocity in the x direction, as in the previous section. This leaves us with a two 
parameter set of solutions, parameterized by e and C,. 



6.3 Perturbative Solutions 

In this subsection we will generalize the work out in Q to the hairy black branes of our system, as 
a function of e and C, at small values of those parameters. Our starting point is Herzog's observation 
that, at e = CXI, the linearized equations of motion about the Reissner Nordstrom black brane at 
1^1=2 admit a regular static solution scalar solution proportional to yq^- explained in 

this solution can be taken to be the starting point for a perturbative expansion of hairy black brane 
solutions in a power series in e. The solutions of |^ were further generalized to nonzero C- 

In this subsection we generalize Herzog's solutions away from the infinite charge limit, to first order 
in a power series expansion in C'(^), i.e to first order in deviations away from the probe approximation. 
This generalization will prove crucial for generating the equations of supcrfluid dynamics including 
effects of back reaction of the supcrfluid on the normal fluid. 

The techniques for obtaining this perturbative expansion are standard. We do not pause to explain 
our computations in detail; in the rest of this section we simply present the results of our calculations. 
As a function of e and ^ (with both taken to be small) we find that the scalar field is given by 



e (21og(r)-log(r2 + l)) 



+ 1 



4r2 + 4 



-2 (r2 + 1) log(r) + (r2 + l) lo g (r^ + l) - 2 ^ ^ 



48 (r2 + ly 



(6.8) 



The functions in the gauge field in (5.2) are given by 



H{r) = (Hoir) + H,{r)e^ + H2{r)e' + ©(e^)) + 0(l/e'), 
L{r) = (Lo(r) + L,{ry + L2(r)e^ + 0(e^)) + 0{l/e') 



(6.9) 
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where 



C(l-log(2)) 



i/i (r) = ^ + ^ ^ lOr-* + 72r4 log(r) - 27r'* log(2) + ISr^ + ISr^ log(2) 

V48(r2 + 1)2 288 (r2 - 1) (r2 + 1)^ 

- 36r* log (r^ + l) - 28 + 45 log(2)^ ^ + O (C'^) , 

iJ2(r) = 7 3^ ^ f 253r^ + 589/ - 589r^ + 48 (r^ + l)^ log(64) 

55296 ((r2 - 1) (r2 + V 

- 336 (r2 - 1) (r^ + l)' log(2) + 576 (r^ + r^) log f;:^^) - 253^) + O , 



and 



(6.10) 



i2M = C'(C). 

The functions in the metric in (|5.2|) are given by 



g(r) = 1 + 1 (5o(r) + SiMe^ + g^{r)e^ + ©(e^)) + (^) , 

jM =0+1 U(r) + ji(r)62 + + O (^1) , 

fcW = + ^ (^oW + fciWe^ + A;2(r)e4 + 0(6^)) + O (^1^ . 



where 



4 (r2 - 1) 2 (r2 ~ l) C , (3r2 + r2(- log(16)) - 3 + log(16)) 
/ol'^j 3r4 3r4 + 12r4 + ^ ' 



-Ir^ + 12r2 - 5 



36r4 (r2 + 1) 

/54r6 + r4(541og(2) - 23) - 36r2(2 + log(2)) + 41 - 901og(2) 

432r2 (r2 + 1) 

+ 18 (3r6 + 3r4 - 9r2 - l) (2 log(r) - log (r2 + l)) j + O (C^) , 
2 (r*^ + - 2r2) (2 log(r) - log {r^ + l)) 



48r2 V 3 (r2 + 1) 

^- 3 ( 576ri" + 989r* + 624r® log(2) - 1538r*^ + 1248r^ log(2) - 1044r'' 

864r2 (r2 + 1) V 

914r2 - 1248r2log(2) + 103 - 6241og(2) 



(6.12) 



(6.13) 
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5o(r) = O (C^) , 
1 



51 W 



6(r2 + l)2 



C (54 (r^ + 1)^-6 (r2 + 1)^ (-9r4 - ISr^ + 3) (21og(r) - log (r^ + l))) (6.14) 

864 (r2 + 1)-* ^ ' 

_ -6r^ - 2lr'^ - lAr^ - 6(r2 + l)^ {r^ + 2r^) (2 log(r) - log (r^ + l)) +4 ,^ 
^ 864(r2 + l)4 + O (C ) , 

jo(r) = 0(C^), 

fr2 - lU (6.15) 



and 

fco(r)=0(C^), 

Mr) ^ rHH-2{r^ + l)lo,irH{r^^ + l)l0E{r^ + !)-!) ^ ^ ^^.^ 
k2ir)^0{e), 

Upon setting ^ = 0, our result exactly matches with the equations 2.30, 2.31, 2.32 in [^,if we 
replace u — l/r in those equations. 

6.4 The first correction to the phase transition curve 

As we have seen above, Herzog's model undergoes a superfluidity nucleation phase transition at | | = 2. 
It is easy to work out how this phase transition curve is corrected at ©(j^). The phase transition 
curve is determined by finding the black brane solution that admits a static regular and normalizable 
solution to the linearized scalar field equations about a Reissner Nordstrom black brane. Now the black 
brane solution, at fixed temperature and chemical potential, depends on the parameter e. In the probe 
limit the metric of the charged black brane is simply independent of the chemical potential and equal 
to that of the uncharged black brane of the same temperature. This metric (and the accompanying 
gauge field) are slightly deformed at C'(^). The condition for the existence of a normalizable and 
regular zero mode to the scalar equation is also, consequently, deformed. 
Let the restriction of the gauge field at infinity be given by 

where Tc is the location of the horizon. It turns out that a normalizable and regular zero mode to the 
linearized scalar equation occurs at 

The temperature of this solution (taking the back reaction of the gauge field into account) is given by 



rp ^ 

TT 
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It follows that ^ at the superfluidity phase transition is given by 



T 



2+3^(2l0g(2)-l) + O(- 



6.5 Boundary Thermodynamics 

Using the solution obtained in the previous section we evaluate the boundary stress tensor charge 
current. For this purpose we use the standard AdS/CFT formulas p8[ 

Boundary stress tensor = T^^ = — i— lim (2 {S'-'Kc.fi-t"^ - K^) - 6(5f; + 

1 



Boundary charge current = = 
Entropy density = s = 
Temperature = T = 



l6nG r->oo 

4G ' 
/'(I) 



lim r^F'^'^ 



(6.17) 



47r.g(l)' 



where and -f^a^ are respectively the induced metric and extrinsic curvature of a constant r surface. 
The result for the stress tensor and current can be parameterized as the form 



rpfll/ 



1 



1 



Au^'u'' + Bn^'n'' + C {n'^u" + ut^n") + 



A-B 



(6.18) 



167rG 

A, B and C are given by the following expressions. 



4 + 2C2 + C*(log(2)--) +0(C« 



59 31og(2) 
144 8 



6241og(2) -451 
13824 



O 



B 



3 3^ 



L + C^ /131 _ log(2) 
36 ^ V432 8 
6241og(2) -451 
41472 



log(2) 



(6.19) 



0(e« 



O 



+ 0(6^) 



04 



While Qi and (52 are given by the following expressions 
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4 + C^ + ic^log(2)-i) + 0(C^) 



131og(2) 493 \ , 
576 27648 / ' 



24 



7 51og(2) 



36 16 



(6.20) 



Using the above expressions we can compute the coefficients in (|2.6|). We find 



167rG {pn) = iri + % [4 + 2^ + 0{Q^)\ + 



167rG(p.) = ^1 [0(C^)] + [1 + 0(C^)] + 0(e4)| + O (^1 



167rG'(P) = r 



4 I ' c 



4 + ^ O(C^) 



7 

36 



(6.21) 



167rG((Z„) = [4 + C'+0(C 



0(6^ 



167rG(<z,) = % [0(C')]+e 



0(e4 



Further the chemical potential and of our solution is given by 

>2 



Moreover we find 



2 M4 4 



__+^2 ^31og(2) . ^ , ^^^^4 



1 

48 



V 32 - 144j+^(^') 
/ _ 71og(2) X X 

I 55296 1152 / ' ^ ' 



2+^ + ^ (-64 



13 , log(2) 



31og(2) 43 



48 V 32 1152 
253 71og(2) 



-O(C^) 



55296 1152 



4G 



log(4)-l^.^^(^4) 



32 



O 



o 



(6.22) 



(6.23) 
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y _ ''c 


+ 




TT 




47re2 




'1 

9 


+ C^( 




[( 


91 




20736 



3 3 M2 3 



log(2) 
108 



(6.24) 



.21) we have verified all the relations 



Using these expressions and the quantities obtained in 
( p^ ) to the order to which we have evaluated our solution. 

6.6 Stability 

In Appendix O we have demonstrated that the solutions we have computed above are unstable when- 
ever I > J A / 1 (working at leading order in ^ and first nontrivial order in ( and e with both quantities 
taken to be of the same order). 

It follows that superfluid flows with C ^ e are unstable. For this reason in the rest of this paper 
we will specialize to the case that C < 0(e). 

It is of some interest to characterize the superfluid phase on a phase diagram with axes labeled by 
y and C (by conformal invariance the phase diagram can only depend on the ratio y). By restricting 
( |6.22 ) to lowest order, and noting that e > 0, we conclude that the superfluid phase exists only when 



1^1 

On the other hand this phase is unstable when 



2 > 



.25) 



'r' - 14 



.26) 



As ^ > it follows that as we increase C at fixed |^|, the superfluid phase first goes unstable and 
then stops existing. It would be interesting to investigate whether this qualitative behaviour persists 
at finite values of 



7. Superfluid dynamics to flrst order in the derivative expansion 

In the previous section we have determined the equilibrium solutions for hairy black branes, pertur- 
batively in e and the superfluid velocity, and separately in an expansion in In this section we use 
the results of the previous subsection as an input into the fluid gravity map. 

The basic idea here is a simple generalization of the ideas spelt out in [^8|, |4^, ^ ^ , 44 , 45[ . 



We search for gravitational solutions that tube wise approximate the 8 parameter hairy black brane 
solutions described in the previous section, with values of the temperature, the chemical potential, 
and m'^ varying in space and time. The tubes in question run along null ingoing geodesies, and 
foliate our spacetime. Technically, this programme is implemented by working in ingoing Eddington 
Finklestein coordinates (as we have been through this paper) but promoting the parameters of our 
solutions to fields that vary in spacetime. 

The fluid gravity map generates the gravitational solutions dual to fluid flows perturbatively in a 
boundary derivative expansion. The zero order ansatz for such a solution is simply the solution (^) 
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with e, Tc, C** promoted to arbitrary slowly varying functions of spacetime. This ansatz of 



course solves the equations of motion (6^), (3^) and (6^), when all parameters are constant, but 
does not solve these equations when these parameter vary in spacetime. As in ^ 43 this 

ansatz may be corrected to obtain a true solution ( systematically in a derivative expansion) provided 
the eight fluid fields that parameterize our ansatz obey certain constraint equations. These constraint 
equations are simply the fluid equations ( [2.5| ) with holographically generated constitutive relations 
for the stress tensor, the charge current, and a holographically generated correction to the Josephson 
equation. 

In this section we implement this programme to first order in the derivative expansion. 
7.1 The method 

As we have explained above, we will search for gravity solutions that tube wise approximate the equi- 
librium solutions of the previous section. In principle our solutions could be labeled by a temperature 
and a chemical potential field in addition to the normal and superfiuid velocity fields. However, for 
calculation purposes we will find it convenient to trade chemical potential for e(a;), the local expecta- 
tion value of the operator O, and a temperature like variable rc{x), together with u'^(x) and C,^{x). 
The precise definitions of our field variables is given by the equations 



p3 



- rl{x)C,{x) 



where (l>{r,x) is the slowly varying bulk scalar field and T^^'^{x) is the boundary stress tensor. The 
functions Ps{x), Pn{x) and ^{x) are given in terms of rc{x), t{x) and C,{x) determined by thermody- 
namics (i.e. from previous sections). As usual, P^"^ — rj^^'^ + u^^u^ . Note that the second equation in 



(7.1) would have agreed precisely with (p.7[) were we to replace the function ^{x) in (7.1) with m^^^ 



As it stands the velocity field defined by (fTTl) is close to but distinct from the velocity field defined 



by dgj) 



The fluid gravity map is generated by solving Einstein's equations tube wise, point by point on 
the boundary. At any given boundary point we can always boost and rotate coordinates so that 

- (-1,0,0,0), rv = (0,l,0,0) 

In the neighborhood of our special point, however, 

= 7«(-l:/3i,/32,/33), = 7n(-«u,l,»^2,?^3), = A + "2/32 + "3/?3 

_ 1 _ 1 (7.2) 

where Pi and rii are of flrst or higher order in derivatives of fluid fields at the special point. 

In this paper we will work only to first order in the derivative expansion. At this order we are 
sensitive only to first derivatives of /3i /32, /Js, n2 and ns along with the first derivatives of ^, Tc and e. 

The solution at our special point preserves an SO{2) symmetry (of rotations in a plane perpen- 
dicular to and n^; the yz plane in our coordinates). This symmetry will help us organize our 
calculation. To start with it will prove useful to organize first derivative 'fluid data', i.e. all the first 
derivatives of the fluid fields at our special point, in terms of their 5*0(2) transformation properties. 
We list our results 
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• First order derivative excitations with spin (scalars): 
^1 = iSiC, S2 = i^ie, ^3 = dol3i, Si = dil3i, S5 = 

Ss = dorc, S9 = dire, S'lo = ^doC, = ^Oqc, S12 = e^jdiPj 

• First order derivative excitations with spin ±1 (vectors): 

^ d,Pi + dil3, [V6],=d,re, [Vj]^^don„ [Vs], ^ d,Pi - di^i 

• First order derivative excitations with spin ±2 (traceless symmetric tensors): 

= dil3j + djPi - {dkPkjSij, [T2]ij = diUj + djm - {dknk)6ij 

Here {z,j} = {2,3}. 

Following the methods of Q , in order to derive the metric dual to a fluid flow we need 

to solve the equations of motion, order by order, in the derivative expansion. That is we set the 
metric g of our solution to gQ + egi . . . (and similarly for the gauge fields and the scalars) and solve 
the bulk equations of motion at first order in e. As explained in |3^, the resulting equations 

are of two sorts. The Einstein and Maxwell constraint equations reduce simply to the equations of 
energy momentum and current conservation, and do not involve the unknown fields gi etc. These 
equations relate some of the independent derivatives listed above to others. On the other hand the 
dynamical Einstein and Maxwell equations allow you us compute the unknown fields gi etc in terms 
of the constrained derivative data listed above. 

7.2 The constraint equations 

We will now first describe the solution of the constraint equations, before turning to the dynamical 
equations. 

In addition to the conservation equations described above, there is one additional source of con- 



straints on the derivative data given in §7.1. Our demand that our solution be asymptotically AdS 
requires, in particular that the boundary field strength vanishes, implying d^^^i, — d^^i, vanishes. We 
must add this equation to the list of equations that constrain independent data. 

It is convenient to decompose the constraint equations according to the its quantum numbers 
under the preserved 5*0(2). We now perform the relevant decompositions, and state which pieces of 
data we use these constraints to solve for. 

• Current conservation: It is a spin-0 constraint. Using this we shall solve for Sn. 

• Stress-tensor conservation: It is effectively four equations. Among them two are spin-0 con- 
straints and one spin-1 constraint. Using this we shall solve for Ss, Sg and Vq,. 

• Curl-free condition on ^^: This imposes a set of 6 equations. 

Two of them transform in spin-0 ([9o'?i — c'l'fo] and Sijdi^j). Using these we solve for Siq and 
respectively. 

Four of them transform in two separate spin-1 ([9i^i — di^i] and [di(,o — 9oCi]). Using these we 
solve for V7 and Vs- 
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After solving for dependent data , the remaining independent one derivative pieces of data are 
given as follows. We have seven spin-0 {Si, • • • , S7), five spin-1 (Vi, • • • , V5) and two spin-2 {Ti,T2) 
boundary data. 

For later use we will find it useful to list covariant expressions for the independent data. These 
expressions are most usefully written in terms of the projector normal to the velocity/ superfluid 
velocity frame 

Using this projector one can write the following covariant expressions for our choices of independent 
boundary data as follows: 



Spin-0 



5i = i(n^9^)C, S2 = ^{n^d^,)t, S:^ = u^n^d^Uu, Si = n^n" d^u^,, 
S5 = P^^a^n., ^6 = P'^'d^u,, S7 = e^^P^n^u^dpU, 



Spin-1 



Spin-2 



[TiU = p;ps [9, 

[T2Vu^P?,P£[d.i 



dpUa 



dpHcr] - SsPfj.^ 



(7.5) 



7.3 The dynamical equations 

Following earlier work on the fluid gravity correspondence ^ , we work in the gravitational 

gauge grr = and grp — Up. For the U{1) field we continue to demand that the scalar field be real. 
With derivatives taken into account this requirement no longer sets to zero, but allows us to 
determine rather simply, by demanding the consistency of the equations for and </>* . 

We will now solve for the first derivative corrections about the basic fluid gravity ansatz. As we 
have determined the equilibrium solutions, in the previous section, only to order we can of course 
compute the metric dual to fluid flows only at the same order in . 

We now describe in rough terms how we determine the deviations away from the zero order fluid 
ansatz. Let us start with the gauge field and scalars. At leading order in ^ we take derivative 
corrections to the gauge field and the scalar field to have the form {6A^ is the derivative correction 



^*As we have indicated above, we solve for some first derivatives of fluid fields in terms of other derivatives. The 
relevant equations are linear and easy to solve; the solutions are explicit but lengthy and we do not present them here. 
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for the gauge field) 



i=l 



o 



(7.6) 



5(j) 



1 



i=l 



o 



We now describe in structural terms how we have solved for these functions, emphasizing boundary 
conditions 

1. It turns out that 5Ai{r) obeys a first order differential equation in r. The general solution of 
SA^ diverges linearly in r while expanded around r = oo. We fix the constant of integration 
(coefficient of the homogeneous solution in this equation) by setting the coefhcicnt of the linear 
term in r to zero. This choice of boundary conditions is forced on us by the requirement that 
the bulk current goes to zero at the boundary so that the boundary current is really conserved. 

2. 5Hi{r) obeys a second order differential equation (arising from the r component of the Maxwell 
equation). The two integration constants for this equation are fixed as follows. One of the 
integration constant is determined from the requirement of regularity at the horizon. The other 
integration constant is obtained from the requirement that there exist a regular scalar field 
solution (see below). 

3. 5Li{r) obeys a second order differential equation given by the x component of the Maxwell 
equation. Here one of the integration constant is determined imposing the regularity of the 
solution at the horizon. The other integration constant is fixed using the fact that according to 



equation (7.1) does not receive any derivative correction. A generic solution of 5Li{r) dies of 



at infinity like -pi; the coefficient of this ^ must be set to zero. 

4. The equation for Xi{r) comes from the y or z component of the Maxwell equation. This is also a 
second order differential equation and its integration constants are determined in a similar way 
as in 5Li{r). 

5. The equation for the scalar field determines 5(j)i{r). Normalizability and the definition of e(x) as 



given in equation (7^) fixes the two integration constants here. More specifically, an expansion 
about infinity of a generic solution to the scalar field equation takes the form 



^0,(r) 



ln(r) 



O 



Our boundary conditions are that both a.; vanishes (from the requirement of normalizability) 
and that bi vanishes (from our definition of e). These two requirements completely fix the scalar 
fluctuation. As described above, the further requirement that the scalar fluctuation be regular 
at the horizon yields a boundary condition on SHi{r) (see above). 
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Let us now turn to the metric field. In the strict limit of ^ — > the scalar and gauge field do not 
back react on the metric. The derivative expansion of the metric in this limit is thus that of uncharged 
fluid dynamics and was determined in Is^] to be 



2uadx^dr + 



dxf'dx" 



1 T UnU, 



P. 



dx^'dx" 



1 



2ru^ {u.d) - - {d.u) + r^F — ct^^ 



(7.7) 



Where 



and 



F{r) ^ -— [-log (r^ + l) +41og(r) - 21og(r + 1) + 2tan"^(r) - tt] 



pa ( dgUp + dpUp, d.u 

The new results of this paper are for the derivative correction to the metric at order O (^). We 
parameterize the corrections to the metric as 



Sids^) 



E - (- 



Ufj^dx^dr 



rdx'^dx'- 



J2 ^' ^f' 

.1=1 



Ui_,Ui, 



-Idx^^dx^l E 



i=l 
5 



E S^ SK, 

{ui^^Uy + n^Ufj) 



(7.8) 



We now describe, very qualitatively, how we have solved for these functions. 

1. 5gi{r), 5fi{r) and dKi{r) are determined solving three coupled equations obtained from the 
(rr) (rv) and [xx) component of the Einstein equations. Once decoupled using appropriate 
combination of these functions, two of the equations become first order and the third one is a 
second order differential equation. Two of the four integration constants are determined using the 
asymptotic AdS condition of the metric. A third integration constant is fixed by demanding 
the regularity of the function 5Ki(r) at r = 7'c. The last integration constant is fixed to ensure 
that the vv component of the boundary stress tensor receives no derivative corrections so that 



the equation (|7.l| ) is satisfied 



^ After this condition is imposed ^ expansion of the functions Sgi{r) and 5Ki(r) take the form 

hm %(r) = O , hm SK,{r) = O 
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2. The function dJi{r) is determined using the (rx) or (vx) component of the Einstein equation. 
This is a second order differential equation in r. The two integration constants are determined 
using the normahzabihty and the definition of boundary stress tensor (according to (^^) the 
vx component of the boundary stress tensor should not receive any derivative correction). The 
general solution for SJi{r) has the following expansion around r = oo. 




Our boundary condition is that jo and ji both vanish. 



3. Yi{r) is determined from the (vy) or (vz) component of the Einstein equation . This is a 
second order differential equation in r. The two integration constants are determined using 
normahzabihty of the metric and the definition of boundary stress tensor (the (vy) or (vz) 
component of the stress tensor should not receive any derivative corrections) . This condition is 
exactly same as that of 5Ji{r) in terms of the coefficients of ^ expansion. 

4. Wi{r) is determined from the {xy) or [xz) component of the Einstein equation . This is a 
second order differential equation in r. The two integration constants are determined using 
normahzabihty and regularity of the metric respectively. A generic solution of Wi behaves like 

at large r. Our boundary conditions are that the leading coefficient of this leading piece 
vanish. 



5. Zi(r) is determined from the {yz) component of the Einstein equation . This is a second order 
differential equation in r. The two integration constants are determined exactly the same way 
as for Wi{r). 



7.4 Results for Bulk Fields 

Without further ado in this subsection we simply present our final results for all the fields defined in 
the previous subsection. 

We have performed all our computations in this section using Mathematica. In several instances 
we have carried out calculations to higher order, in the Mathematica file, than we have presented 
below, mainly to avoid burdening the reader with very lengthy expressions. 

The solutions presented in this subsection determine the full first order correction to the gauge 
field, scalar field and metric to the relevant order in an expansion in e and ^ . Now we choose to scale 
C like e. We present our results below in terms of the order one field 

c 

X = - 

e 



Recall that, according to the results of §3.6, our fiuid becomes unstable whenever x exceeds a number 



of order unity. So while x can be arbitrarily small, it is unphysical for x to be made arbitrarily large. 
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Results for the gauge field and scalar field 



Mi(r) = e 
M2(r) =e 
6A3{r) = - 
SA^r) = 
SA^ir) = e 



(96x2 - 5) + 48x2 + 1 



14^3 



(2-3r2)x 

(2r2 + 1) X 

16 (2r2 + l) x^ 2r 
7^ 

(1 - 5r2) X 



14^3 



3(r2 + l) 



3 (r2 + 1) 



7r3 



(7.9) 



SHi{r) 
SH2{r) 
5Hs{r) 
SHiir) 
SH^ir) 
SHe{r) 



r(r + 2) (96x^-5) -48x^-1 



= e 



14r(r + 1) (r2 H 

(3r2 + 6r + 2) X 
7r {r^ +r^ + r + 1) 

(5r2 + lOr + l) X 



1) 



+ 0(e3) 



7r (r^ + + r + 1) 
16 (2r2 + 4r - 1) x^ 



7r (r^ 



= e 



(5r 



-r + 1) 
10r+l)x 



14r (r^ + r2 + r + 1) 
[2r(r + 2) - 1] x^ 



+ O(e') 
0(e3) 



7r(r + 1) (r2 + 1) 



(7.10) 
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X (log + l) - 21og(r + 1) + 2tan-i(r) - tt) 



4^2 



+ 0{e 



= e 



log (r^ + l) - 21og(r + 1) + 2tan-i(r) - tt 



+ C(e 



0(e 



X (log (r^ + 1) - 41og(r) + 21og(r + 1) - 2tan-i(r) + tt) 



3r2 



+ 0(6^) 



0(e 



X (log (r^ + 1) - 41og(r) + 21og(r + 1) - 2tan-i(r) + tt) 



6r2 



(7.11) 



X2ir) = 



log (r^ + l) - 2 log(r + 1) + 2 tan"^ (r) - n 



X (log 



96r2 

+ 2tan-i(r) - tt) 



4^2 



+ 0(e* 



X3(r)=0(e4) 



(7.12) 



^5(r) = e 



X (log (r^ + 1) - 4 log(r) + 2 log(r + 1) - 2 tan-i(r) + tt) 

4^2 
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A(i-8x2) (tan-i(r-)-log(l + r)-^) 



+ - (36x' - 1) log(r) + ( - - 6xM log (r^ + l) 



-— X (-7 log (r^ + 1) + 41og(r) + 101og(r + 1) - lOtaii-i(r) + Stt) 



i(r) = e 



X (-7 log (r2 + 1) + 8 log(r) + 6 log(r + 1) - 6tan-i(r) + Stt) 

14 (r2 + 1) 

4x^ (-7 log (r^ + l) + 121og(r) +21og(r + 1) - 2tan-i(r) + tt) 



7 (r2 + 1) 



X 



i>5(r) = -4 (-"^log (r^ + l) + 81og(r) + 61og(r + 1) - 6taii-i(r) + Btt) 



Results for the metric 



-X^ (-71og (r^ + 1) + 121og(r) +21og(r + 1) - 2tan-\r) + tt) 



(7.13) 



'5/i(r) 




"-2(80x^-3) 


= e 




<5/2(r) 


= e 


16x" 
21r4_ 


+ 0(e3) 


Sfsir) 




ri2x 

' 7r4 


+ 0(e^ 



Sf4{r) 



+ 



320x2 _ (r^ + 1) (-")los (/■- + 1) + 8los(/-J + 21os(/- + Ij + 4t,au- ^(r)) 
' 21r4 

2 (r^ + 1) (r^ (7rr2 - r + tt - 3) - 2) 

6 _L 



9r2 



9 (r6 + r4) 



(7.14) 



+ 0(e2) 



<^/5(r) = e 

'5/6(r) = - 



6x 



7^4 



+ 0(e3) 

320x2 (r* + 1) (-5 log (r^ + l) + 8 log(r) + 2 log(r + 1) + 4 tan"! (r)) 



21r4 



+ 



2 (r^ + 1) (r2 (7rr2 - r + tt - 3) - 2) 
9 (r^ + r^) 



9r2 
+ 0(e2) 
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<55i(r) = 0{e^) 



5gi{r) = 



18 



(-5 log (r^ + 1) + 8 log(r) + 2 log(r + 1) + 4 tan'^ (r) - 27r) 



+ ■ 



9(r + 1) (r3 +r)^ 



(7.15) 



^56 (r) 



1 

~2 

+ ■ 



18 



(-5 log (r^ + l) + 8 log(r) + 2 log(r + 1) + 4 tan'^ (r) - 2it) 



J.6 + 4^5 _^ 4^4 _^ g^3 4. ^2 _ 2^ _ 2' 
9{r + l){r^+rf 



5K2{r) = 0{e^) 



— (-Slog (r^ + 1) +81og(r) +21og(r + l) + 4taii-^(r)) 
o 



+ 



4 - 2r2 (Trr^ 



r + TT — 3) 



3 (r2 + 1) 



+ 0(e') 



(7.16) 



5Ke{r) 



= Oie^) 


1 


V2 


~ ~2 


3 




4 


+ — 



(-Slog (r^ + 1) +81og(r) +21og(r + 1) +4tan-^(r)) 



4-2r^ Trr^ 



TT - 3) ■ 



3(r2 + l) 



+ 0(e2) 
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SMr) = -e' 



(r^ - 1 X 



^^32^ (log (r^ + 1) - 21og(r + 1) - 2tan-i(r) + tt) 



+ ^(2-3r) 



1) (-log (r^ + l) +21og(r + l)-16tan-i(r)+87r) 



96r2 



+ ■ 



27^4 _ 3^.3 + 20r2 - 3r - 19 



5J3(r) = -e^ 

W4(r) = 0(e3) 
5J5(r) = 0(e3) 
W6(r) = 0(e3) 



6 (r2 + 1)' 



144 (r^ + r) 



Yi{r) 



12 (r) = 



+ 



96 r2 

27^4 _ _^ 20r2 - 3r - 19 



[-log (r^ + 1) +21og(r + 1) - 16tan-i(r) + Stt] 
+ 0(e^) 



144r(r2 + 1) 



— e 



(r" - 1) X 
4r2 

6 r 



(log (r^ + l) - 2 log(r + 1) - 2 tan'^ (r) + tt] 



(2 - 3r) 



r3(r) = 0(e^) 



n(r) = 



6(r2 + l)2 
l5(r) = a(e^) 



W2{r) = 



-3iT (r^ + r) + 6 (r^ + r) tan-i(r) + Gr^ + 4' 



16 (r3 + r) 

Stt (r^ + r) + 6 (r^ + r) tan"i(r) + Gr^ + 4' 



W3(r) = 
W4(r) = 
WsW = -- 



32 (r3 + r) 

X [-Stt (r^ + r) + 6 (r^ + r) tan-^(r) + Qr^ + 4] 
32 (r3 + r) 

X [-Stt (r^ + r) + 6 (r' + r) tan-i(r) + Sr^ + 4] 
32 (r3 + r) 



+ (!^(e5) 
+ 0(e^) 

+ 0(e4) 



- - ^ + 5 log + 1) - 8 log(^) - 2 log(r + 1) 

-4tan-^(r) + 27r + O(e^) 



(7.17) 



(7.18) 



(7.19) 
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Zi{r) 




+ ^ - 5 log (r^ + 1) + 8 log(r) 



Z2{r) 



+ 21og(r + l)+4tan-i(r)-27r +0{e'^) 

3 r-37r (r^ + r) + 6 (r^ + r) tan-^(r) + Sr^ + 4 



(7.20) 



16 (r3 + r) 



7.5 The C limit 

The gravitational solutions presented above are complicated largely because they possess very little 
rotational symmetry. At any given spacetime point we have a normal fluid velocity and an independent 
superfluid velocity. These two velocities together break the local Lorentz group at a point down to 
the abelian group SO{2). While we have usefully organized the results of our gravitational calculation 
in representations of 50(2), as representations of S0{2) are all one dimensional, our solutions admit 
several different functions of r. 

In the special case that C = 0, however, the residual symmetry group about a point is SO{2i). 
SO{3) representation theory is considerably more constraining than SO{2) representation theory. This 
implies that the gravitational dual to superfluid dynamics should be considerably simpler in the special 
limit C ^ than in the generic case. 

Let us first present a brief ab initio analysis of the nature of the gravitational solution when C = 0- 
All independent first derivative data may be organized into SO{3) scalars, vectors and tensors. These 
may be chosen as follows 



Note of course that an 50(3) vector or an 50(3) may be decomposed into an 50(2) vector and 
a scalar, while an 50(3) tensor is composed of an 50(2) tensor, vector and scalar. In 50(2) terms, 
therefore, the data listed above totals to 7 scalars, 5 vectors and two tensor. 

It follows from symmetry considerations (and the fact that our parity conserving gravitational 
system will never generate a parity violating vector term, so we can ignore the third vector above) 



Scalar 



d^u^ and P'^'S^C 



Vector 



u^^d^u", P^^a^e and e^""^" u^dxCa 



Tensor 
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that it must be possible to write the metric and gauge field, in the C ~^ limit, in the form 



ds^ = -2.9 u^'dz'^dr + 

+ rlF (—] an^dx^'dx" 



dxf'dx" 



Up,dx^dr 



u^Uijdx^dx'^ 



u^dx^dx'^ 



Ti 



.(c) 



dx^'dx" \ + O 



(7.21) 



A^—H 



Ax 



dr 



+ ■ 



{P'"'d,C^) 



{u.d)u^df, + ^C2 ' 



P'^'d.eda + O 



The results of the previous subsection must obey several relations in the limit C — )■ for them 
to agree with the form presented in ( 7.21 ). In Appendix |d| we have explicitly verified that each 
required relation is indeed obeyed. The results presented in the previous subsection are consistent 
with the form (7.21) once we make the identifications 



20A direct comparison between these two forms is complicated by an irritating feature; the coordinate choice of the 
irevious subsection differs from the one above (it breaks manifest 50(3) invariance) even in the limit — > 0. In Appendix 
J we have explicitly performed the coordinate change that allows one to transform the results between coordinates. 
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Vi(r) 



3(r2 + 1)2 



V2(r) 
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2(r + l) 4 
+ 1 


+ 5 log (r^ 




- 4tan~^(r) + 27r 
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T2{r) 
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= 0, (Required by Weyl 


invariance 










^2(r) 


= Oie"^ 
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■^2(r) 


6e 









;log(r) - 21og(r + 1) 
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and 



- 14.3 
Hi(r) = 0(62) 

^2(r) 



2r 



3(r2 + 1) 



+ 0(e2) 



-5r(r + 2) - 1 



0(e3) 



(7.23) 



Ur{r + 1) (r2 + 1) 

log (r^ + 1) - 21og(r + 1) + 2tan-i(r) - tt 



>C2(r) 



96r2 



-0(e4) 



7.6 Stress Tensor, Charge current and the Josephson Equation 

The results of the previous subsection may be used to read off the values of the boundary stress tensor, 
the boundary current and the correction to the Josephson equation at first order in the derivative 
expansion. Like all the calculations in this paper our results are obtained in a power series expansion 
in e and 

We parameterize our boundary stress tensor and current as 



1 



167rG 



1 



167rG 



(7.24) 



where A, B, C, Qi and Q2 are functions of e{x), C(a;) and rdx) as given in equations (6.19) and 



(3.2c). We further expand the corrections to the perfect fluid stress tensor and current as 
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(7.25) 



i=l ^ ^ 



Our results are given as follows. 
Results for stress tensor: 

ii = 0{e^), t2 = 0{e^) 
Results for current: 



ai = e 



^(3-80 x' 



8x 



1 

24 



03 



04 = 



160 



+ 0(e2), 64 = 0(e^) 



9x 



80 



(7.26) 



(7.27) 



ci = ^ + C9(.^) 
C2 = e2x + 0(e^) 
C3 = O(e') 
C4 = O(e') 

c. = e3x(^i±^KC)(e^) 
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Results for the correction to the Josephson equation: 



+ 0(e3) 



St^3 



32 x' 
7 



(7.28) 



5 % 



/16 X' 



V 7 



In C ^ limit this derivative corrections to stress tensor, charge current and the phase equation 
take the following form 

Oie') 



lim ^f"' = -2rfa'"' 
lim fidii 



(7.29) 



where 



ai 



7 ' " " 24 
7.7 Weyl Covariance of our bulk fields and boundary currents 

In this section we will demonstrate that our fluid dynamical solutions must, on general grounds, obey 
certain constraints that follow from the requirement of Weyl invariance. We then verify that our 
explicit solution does indeed obey these constraints, providing a nontrivial check on these solutions. 

All computations reported in this paper have been performed for superfluid motion on a flat 
boundary metric. However our final results must be the restriction to a flat boundary of results that 
apply in a general weakly curved space. The (boundary) generally covariant version of our final bulk 
metric, stress tensor etc are all given simply by promoting all derivatives to covariant derivatives 
(ambiguities in this procedure and boundary curvature terms all show up only at second order in the 
derivative expansion). 

Given these results in a general boundary spacetimc, it follows on general grounds (see |^) that 
our bulk metric, gauge field and scalar fields must enjoy invariance under the following spacetime 
dependent Weyl transformations and coordinate redefinitions. 

^■tp{v,Xi} ~ _ ^ p-2i/j(t),a:i) 



c = c, 



Note that the Weyl transformed metric g^^ is, in general, not flat even if the original metric is. Let 
us work in the special case that the original metric g^^ is taken to be flat. The boundary connection 
with respect to g^i, is non zero and is given by 
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The new frame covariant derivatives of and arc given by 

^tj.Uu = e"'^ [d^Ui, + Uf^di,ip - ri^„{u.d)ip] 

Using these expressions one can deduce the transformation properties of the scalar, vector and 
the tensor forms appearing in the bulk solution 



S3 = e'^ [^3 - (n.9)V'] , ^4 = [^4 - (u.d)^] , 
S5 = [S5 - 2(n.a)V] , ^6 = e* [Se - 2(u.a)V] 



Vi 



V3 



V4 



(7.30) 



To 



If wc transform the g;mge field and the metric from the new Weyl frame (the frame with tilde variables) 
to the old Weyl frame (the frame where the variables are denoted without tilde), the equilibrium 
solution itself generates some new terms due to the r coordinate redefinition. In the new frame the 
coordinates are r = re'^i'"^^^) and x^^ = x^. This implies the following transformation rule for the 
differentials. 



dr = e'^'f"'^') {dr + rdx^d^^) 
dx^" = dx" 



dr 

dn = -^dr + df, 
ox^ 



(7.31) 



This induces the following transformations on gauge field 



A= — 



re 

r. 



1 

+ — 



(7.32) 



= -re'^ [H (r) (u.SV') + L (r) {n.di})] dr + r[H (r) {u.d) + L (r) {n.d)] 
In the last line we have used the scaling symmetry to set Tc = 1. 
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Similarly the equilibrium metric also transforms and the nontrivial transformation is generated 
due to the term dx^dr. 



2g [ ^ j Uf^dx^dr = -2g ^yj Uf, dx^" {dr + rd^-ipdx") 
-2g [ — j u^dx^dr — 2rg ( — j u ^Ui,{u.d)^]dx^ dx'^ 



r 

rg I — 



dx^dx" 



(7.33) 



-2g (r) u^dx^dr — 2rg (r) u ^Ui,{u.d)'4)dx^ dx^ 



rg (r) 



dx^'dx" 



Here also in the last step the scaling symmetry is used to set Tc =1 



Combining these transformations we find the transformed metric, gauge field and scalar have the 



expected form (expected according to (7.7) ) together with some additional pieces that multiply a 
single derivative of ip. The coefficients of these unwanted pieces themselves have no derivatives, and 
must vanish in order that our result respect Weyl invariance. This requirement imposes the following 
simple algebraic conditions on the fields in the metric, scalar and gauge field: 



SAsir) + 2 SA^ir) - rL{r) = 

M4(r) + 2 6Ae{r) - rH{r) = 

SH3{r) + 2 SH5{r) = 0, SH^ir) + 2 SHe{r) = 

SLsir) + 2 SLr^ir) = 0, SLiir) + 2 SLe{r) = 

6(b3 (r) + 2 (505 (r) = 0, (504 (r) + 2 SMr) ^ 

Xair) - Xiir) ^ 



(7.34) 



Sh{r) + 2 6h{r)^0 
SUir)+2dfeir) + 2rg{r)=0 

SJsir) + 2 SJ5{r) + rg{r) = 0, SJiir) + 2 SJe{r) = (7.35) 
SKsir) + 2 SK^ir) = 0, SKi{r) + 2 SKeir) = 
Y3{r) ~ Yiir) ~ rg{r) = 0, W^ir) - Wi{r) = 



We also require that the stress tensor, charge current and Josephson equation in our model are 
invariant under Weyl transformations. As these boundary quantities are all independent of r, the 
redefinition of r is irrelevant to the study of Weyl transformations of these quantities. Using only the 



equations ( 7.30 ) we find the following constraints on the coefficients in ( 7.25 ) 



P3 + 2 P5 = n + 2 Pe = 
V3 -Vi =0 



03 + 2 05 = 04 + 2 ag = 
63 + 2 65 = ^4 + 2 fee = 
C3 - C4 = 

6^13 + 2 (5/i5 = (5/i4 + 2 Sjie = 



(7.36) 



The equations ( 7.35 ), ( 7.34 ) and ( 7.36 ) must apply to any consistent asymptotically AdS solution 
of gravitational equations. In particular these equations must apply to the results of this paper, and 
constitute a nontrivial consistency check on our algebra. We have explicitly checked that the results 
of §7.4 and 37.61 obey these constraints, to the calculated order in e and 



7.8 Entropy Current from Gravity 

Fluid flows obtained from the fluid gravity correspondence are automatically equipped with families of 
local entropy currents of positive divergence. A particularly natural choice for this entropy current was 
presented in equation 3.11 of [^9|. Using this formula for our solution we have computed the entropy 
current dual to our fluid flow. This entropy current has a piece at C(l) and a piece at ©(l/e^), and 
takes the form 



4GJr = r,^u- + ^^5.(«r^ 



(") 



(7.37) 



where 



(u 
{u 

Kg 

K4 

(u 



7 ^ 
7 



(n) 



_i /18X 

2 V 7 

1 /160 

^ "2 7 



1 



r 



9X 
7 

1 1^80 



80 x') 
+ 0{e'), 4"^=0(e^) 



e^X + 0{e^) 



(7.38) 



4") = Oie') 



Ki 


_ 1 

~ 2 24 






1 2 

= r ^ 


+ 0(6^) 


(v) 

Kg 


= 0(e3) 




{v) 
K^ 


= 0(^3) 




{v) 
«5 







(7.39) 
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Quite remarkably this gravitational entropy current agrees exactly with the simple fluid dynamical 



current described in ( 2.21 ) (to the order at which we have done the calculation). 

Actually, all but the first two terms on the RHS of (2.21) are 0{e^) or higher. It follows that to 



•^5 



r-p diss 



(7.40) 



We have checked that the gravitational entropy current, presented in this subsection, exactly agrees 
with this form to this order. 



8. Transformation to standard frames and identification of the dissipative 
parameters 

The equations of gravitational super fluid dynamics, derived in the previous section are presented in 
a modified phase frame that is adapted to the expectation value of the operatore e(cc), and is not 
particularly natural from a fluid dynamical point of view. 

In this section we will transform our results to the ^diss — fluid frame and the transverse fluid 
frames. We will then compare these results with the general 'theory' of dissipative dynamics presented 
in We will find perfect agreement with the general structures predicted in and so be able to 
read off the values of all 10 nonzero dissipative fluid parameters. 

To begin this subsection we first recall the structure of the boundary equations that emerged 
from our gravitational computations of the previous section. In the previous section our gravitational 
solutions were parameterized by the eight fields u'^{x), e{x), rdx) and C'^{x). In this section we will 



find it convenient to use the first of (6.22) to define a new field fio{x) and to eliminate e{x) in favour 
of the new field Ho{x). We will always suppose that this has been done in what follows. 

Our gravitational results for the stress tensor, current and superfluid phase are of the form 

T''" = [Pnirc, Mo, ^o) + Pirc, ^^oAo)] W'u'^ + P{r,, ^io, Co)??^" + /(Co, A^o, r-c)Co^Co + ^^"^ 
J'' = g„(r„eo,A^o)u^ - fir,,Co,^io)C^ + J^s (8-1) 

where the functions p„, P etc are the thermodynamical functions derived in §^ Our results are 
presented in a modified phase frame. We remind the reader that this means that we have presented 
our answers in terms of a new auxiliary phase field and its modulus defined as 

= -^i^u^ + r,C^ 



^0 = Va^o-^^C^ 

Note that is not equal to the phase field (because /iq ^ mO- Instead the relationship between 
these two fields is given by 

= Co^ - ^^d^ssU^' 

c _ c /^o (8-3) 

? — ^0 ^ IJ'diss-T- 

In this section we wish to transform our gravitational results into the fidiss — fluid frame and 



the transverse fluid frame. As we have explained above, (B.l) is in a modifled phase frame not a fluid 
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frame. This means th at (|8.l| ) differs from the form (2.6) in two ways. First thermodynamical functions, 
hke the pressure, in (|8.l| ), are functions of {fJ.o,rc,£,o). According to the form specified in (2^) the 
choice and definition of two of the three thermodynamical variables- e.g. the chemical potential and 
temperature - is upto the user, and can reasonably chosen in a fluid frame to be /io and rc- If we work 
in a fluid frame, however, the third variable has to be ^, the magnitude of the phase field. As we have 
seen from ( ^.3[ ) 7^ 

The second way in which (8A) differs from the from the standard fluid frame form (2^) is that 
the fourth term in the expression for T'^'^ in (B.I) is proportional to ^q^o- Similarly the second term 

in 



in the expression for J^- 



diss 



.1) is proportional to ^g. The corresponding terms in (2.6), however. 



are proportional to ^^■^'^ and respectively. 

These discrepancies are easily cured. In order to move to a fluid frame, all one needs to do is 



substitute the expressions for and as functions of ^ and ^'^ (see (8.2)) into (pH]). That is we 
must perform a prescribed field redefinition on ^g to take it to the field In addition we are free to 
perform any additional field redefinitions 



rc = fc + 6rc 



(8.4) 



in order to transform to any particular fluid dynamical frame. Of course (5/i, 5rc and 5u^ above are 
all necessarily of first or higher order in derivatives. 

We will now describe how these general ideas can be used in practice to transform our gravitational 
results into the ^diss = frame and the transverse frames respectively. 



8.1 Transformation to the ^idiss ~ frame 



In order to transform to the ^diss — frame we must take 5/i — —^diss in (3-4). We can immediately 
work out the effect of this field redefinition, combined with the change of variables from ^g to ^'^ , on 
the current and stress tensor. Wc find 



and 



= Qnirc, Co, A^o)"^ - f{rc, ?o, ^0)^0^ + J^s 

= q„(r„e,Ai)^^-/(re,e,M)C'' 

+ {dq.n + A^od/ )u^ + ^id^ssf{rc, ^0, A'o)"^ - df ^ + 



T^" = [Pnirc, Mo, Co) + P(.r,, Mo, Co)] W'u'^ + P(r„ ^o)v'"' + /(Co, Mo, ^c)eo^Co" + 
= [Pnirc, M, + Pirc, M, 0] u^'u'' + P{r,, /i, ^v'^'' + f{r,, 
+ [dp + fil df] u'^u-' + ^Xd^ssf{rc, Co, Mo)(«^eo + W'eo) 
+ df CC + dP P^"" + ^^'^ 



(8.5) 



(8.6) 



where /(r^C, m) 
given by 



gs(rc^^4i) _ pAra^i,tJ.) ^-^^ operation 'd' acting on any function of Tj,, S, and /i is 



dA = -fidi 



Mo 



r9A(re,e,M)1 




raA(re,e,M)l 


[ \ 


f^diss 


9m 
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The gravitational stress tensor and charge current have now been recast into the fluid dynamical form 
with 



TT^- = [dp + pI df] u^^u- + Md.ss/(rc, Co, + u-i'i,) + d/ C^C" + dP P^-^ + tt'^ 

f^diss — 



■7) 



While TT^'^ was orthogonal to the velocity field u'^, the same is not true of tt'"^ in ( p.7[ ). In order to 
enforce the transversality of tt^"^ (a defining condition for the fj.diss = frame), we must now redefine 
Tc and m''. The determination of Svc and (5m'' is particularly simple at leading order in (the order to 
which our gravitational results have been obtained). Recall that, to leading (unit) order in tt^''' is 
already transverse to u^. The piece of T'^" that is not transverse to starts out at O (^). It follows 
that the 6rc and (5m^ must both be of order 0{-i^). It is important to implement this field redefinition 
only in the 0(1) part of the equilibrium or perfect fluid part of T^"^ (as we are keeping track of the 
final answer only 0{-ij) and only to first order in fiuid derivatives). This piece and its transformations 
are given by 

equilibrium at O { j 

=ri{Au^u'' + r^n (8.8) 



=f:i(4{t^{t'' + 77^'')+4f^ 



uT 

— (Sm^u'' + P''") + (iJu^'m" + (5m"u^) 

fc 



Adding (B.8) to ( |8.6| ) gives the complete expression of tt'^" after this variable redefinition, tt^'' is 
transverse if we choose 



Ur^Src = - {dpn + pi df - 2p)f p^i 



(8.9) 



(as can be verified by dotting tt^'' with u^u'^ and with u^n'^). With Svc and chosen as above 

= df Qf'C + [dP + 4rf 5rc) P^" + ^'"'^ (8.10) 

(all expressions involving cancel). 

Note that the corrections to J'^^g^ that arise from the field redefinitions of r^, and are all at 
0{-^). As we have not kept track of the current to this order in our gravitational computation, we 



will ignore all such terms; Jl^^^^ continues to be given by (8.7) 



Note that the final expression ( B.IO ) for tt'^'^ depends on Sr^ but not on 5u^ (this is, of course, true 
only to first order in the derivative expansion). While Stc is given by ( |8.9[ ), it is equally well given by 
solving the equation 



df C + ^ {dP + -ir^Src) = 



.11) 



That (3.11) must be true follows from the trace of (5.1C) and the observation that tt'^" and tt^'^ are 



both traceless. That the expressions for Src obtained from (B.ll) and (8.9) agree follows from acting 
the operator 'd' on the equation 

i-Pn +3P- fe) = 
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which itself is simply an assertion of the tracelessness of the perfect fluid stress tensor 



Inserting the expression for 5rc from (3.11) into the expression for vr^"^ in (8.10) we conclude that 



the gravitational expressions for tt^'' and Ji'j„„ in the /i^iss — frame are given by 



die 



.13) 



In the rest of this subsection we will now compare this gravitational result to the 'standard form' 
for the dissipative corrections to the stress tensor and charge current predicted by our 'theory' of Weyl 
Invariant fluid dynamics in §^ This standard form was given in (3.18) and is reproduced here for the 
convenience of the reader 



[PsS + PbSb + P^S^] n^'n" 



{QsS + QbSb + Q^S^\ 
+ [RsS + RbSb + Ry,S^]n^' 

+ CaV,^ + CbV^ 



(8.14) 



where 



and 



5 = 



Syj — Tl^Tl (7 ny 



The standard form described above uses a basis of independent first derivative terms that differs from 
our choice of independent first derivative terms in our gravitational solution. In order to compare 
the two forms we simply reexpress Sa^ Sb and Va in terms of the independent first derivative forms 



= d{-p„ +3P- ft) = -dpr, + SdP - fdf - 

: C^df + 3dP - [dp„ + i?df - 2pf 



: -dp„ + 3dP - fdf + 25/ 



(8.12) 
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employed in the gravity calculation, utilizing the equations of motion. We get 



5d 



(8.15) 



We then compare the expressions so obtained with our gravitational results ( ^.ISj ). We find that 
the two expressions match perfectly if we choose 



167rGQf, = 
167rGQ^ = ^ [240x' + O(e 



4 



-240- + 0(e)° 



4 



167rGi?5 = 



TT 
72 



-240- 



O 



167rGi?„ = ^ [288 ex^ + 0(e)'] + 



(8.16) 



167rGP6 
167rGPu, 



[240x' + O(e)] +0 

0{eY] +0 



ex 

^3 



-6 



3x' - 288 X 



0(e)^ 
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167rGa = 



-l + log(4) 



^3 r 



167rGT = -27r-^ 



-4 



^3 r 



-4 



) 



\ ^ 



0{ef 



O 



O 



(8.17) 



Note that — Pa, ~ Pt, Ra ~ Qb and Cf, — Ea, so that our resuhs obey the Onsager 
relations hsted in §^.5[. It may also be verified that these results obey all the positivity constraints 



required on general grounds in § 3.5 



8.2 Transformation to the Transverse Frame 

Following the discussion of the previous subsection, it is also possible to cast our gravitational results 



into the transverse fluid frame. As mentioned in §3.3, the basis of first derivative quantities with non 
zero coefficients most suitable for this frame are 



(8.18) 



These quantities may be expressed in term of the quantities (defined in (7.3)) used for the gravity 
calculation as follows 



1 



e- — {2Si^Se)+e 



28 

2S4 — Sq 



14 



28 



2 T^'^ 



(8.19) 



Let us rewrite the first derivative corrections to charge current obtained from gravity (given in ( [7.25 ) 
and ( 7.27 )) in the following schematic form 



1 



IQttG e 



(8.20) 



In the gravity solution the stress tensor ( tt'^") is given as the following (see ( [7.26 )) 

167rG ^^''^ = [-2rln^"' + ©(e^)] + O ' ^ 
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Then using the procedure outhned in the previous section and in §(2.4) we can compute first 
derivative corrections to stress tensor, charge current and chemical potential in the transverse frame 
(which we denote by 7r^"^\ (-^dS)/^ ^'^'^ A'dlss respectively). We find 



where iidiss is the first derivative correction to the chemical potential obtained fi'om gravity given in 



(7.25) and in (7.28). Here the formulas presented in (3.21) are valid only at the leading order in e for 
each independent data at one derivative order. 



As in the previous subsection, we then consider the expected standard form fiuid expression which 
is given by 



{PaSa + PbSb + PwSw) ynf_,n^ - 
+ Ea {V^n- + V^n^) + Eb {V^n^ + V^n^) 



{RaSa + RbSb + RvjSw) 

+ CaV^ + CbV^ 
IJ'diss = [QaSa + QbSb + QwSw] 



.22) 



Just as in the previous section the gravity result after the frame transformation (8.21) perfectly 
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fits into the above form provided we identify equation 



and 



Qa = 16ttG 



Ra 
Pa 
Qb 
Rb 
Pb 

Qw 

Rw 
p 

in 



25e 



52 



+ 0(e)° 



O 



0{ef 



24 
25 



O 



z5e 



(167rG)e2 
(167rG)e2 



288 Q n9 



(167rG)e2 
37r3 



288 Q n9 



IGttG (167rG)e2 



IGttG 
167rG Ca 



[0{e 



O 



WttG Eb = -27r3 



IGttG Gh = 



6 



-l + log(4) 



-2x^ e2 + 0(e) 



0{er 



O 



^3 r 



167rG T = -27r^ + 



6 



2x^6^ + 0(6^ 



.23) 



.24) 



Note that in this transverse frame also we have Q^j = Pa, Rw ~ Pb, Ra = Qb and C'b ^ Ea, which 
constitutes the expected Onsager relations. All the positivity constraints given in §3.4 are also obeyed 
in this frame. 

In — )• limit derivative corrections to stress tensor, charge current and the phase equation in 
transeverse frame take the following form 



lim 

lim 



AT) 



j{T) 

diss 



.25) 



lim iSF^ ^ 

A ^c\' diss 
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where 



/3i - -271^ + ^ 



6 



1 



IQttG 
IQttG 



7r3T3 



52 



.26) 



8.3 Transformation to the Landau-Lifshitz-Clark-Putterman Frame 

Our results may also be transformed to the Landau-Lifshitz-Clark-Putterman frame. We have not 
fully worked through the complicated algebra needed for this process. However we have verified that 
this transformation yields 



167rG g^'^ = ^ 



oil 



Note in particular that is not zero, predicted by the formulas of Clark and Putterman. 



9. Discussion 

Our paper suggests many directions of for future research. In this paper we have studied the bulk dual 
to Einstein-Maxwell systems in the absence of a bulk Chern Simons term . A Chern Simon's term has 
already been shown to lead to yield qualitatively new contributions to charged fluid dynamics even in 
the absence of charged scalar fields |^9[ ^ . It seems certain to contribute new and interesting 

terms to charged superfluid dynamics as well. We leave a study of the effect of this Chern Simons 
term on holographic superfluid flows to future work. 

Relatedly, we should emphasize that the framework for dissipative fluid dynamics presented in this 



paper crucially assumes that the superfluid entropy current takes the canonical form discussed in §2.5 



While this assumption turned out to be true of the gravitational system studied in this paper we do 
not know of any physical reason for it always to be true. In particular we think it very likely that the 
entropy current will be modifled in superfluids dual to a gravitational system with a nonzero Chern 
Simon's term as described in the previous paragraph. It would be very interesting to understand the 
general rules for the entropy current at flrst order in the derivative expansion, and thereby to construct 
the most general framework for flrst order dissipative superfluid dynamics. 

It would of course be interesting to study dissipative corrections to holographic superfluid dynamics 
away from the rather strange limit (of very large e and perturbatively in the superfluid condensate) 
employed in this paper. Infact it is interesting that a continuation of the model studied in this paper 
to a particular order one value of e and supplemented with a particular Chern Simons term and a 
more complicated scalar potential (however with squared mass =-4 as in this paper) is a consistent 
truncation of IIB supergravity on AdS^ x S'^ ||5l|. It would be fascinating to be able to work out 
the map from gravity to fluid dynamics in this system, though that might require numerical work. 
It would also be interesting to follow the investigation of to study whether hairy black holes in 
global AdS'i (see e.g. ^ ^) can be thought of as stationary superfluid flows. 

The task described in the previous paragraphs is hampered by the lack of explicit analytic solutions 
for equilibrium conflgurations. One could imagine circumventing this lack in two different ways. First 
it would be very interesting to attempt an abstract analysis of dissipative fluid flows, along the lines 
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described in Appendix |^ (following Sonner and Withers) for equilibrium flows. The aim of such an 
abstract analysis could be to prove that the gravitational entropy current always agrees with the 
canonical form presented in ^ in the absence of bulk Chern Simons terms, and to prove that all 
Onsager relations are obeyed. This would constitute a proof that gravitational super fluid dynamics 
falls within the framework of dissipative superfluid dynamics described in §^ above (in the absence of 
bulk Chern Simons terms) . A second possible approach to the same problem is numerical; a clever use 
of numerical techniques should permit the numerical generation of the gravitational solutions dual to 
arbitrary first order superfluid fiows at finite values of e. 
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A. More about the canonical entropy current 

A.l Frame invariance of the fluid frame entropy current 

In this subsection we will demonstrate that the entropy current 



is frame invariant. 

Let us perform the following change of variables (we work accurately only to first order in the 
derivatives) 




diss 



T 



T = T + ST 



(A.l) 



ji = ^ + Sfx 



Under this redefinition Jj^^^g and tt^^ change according to the formulae given in 



(Hi- 



d{ps)+v'""dP 



(A.2) 



diss 



gnSuf^ + dqnU^ - dqs — 



In the new frame the entropy current is given by 
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-,- ^.jt. . (tt^'" + fa^") 
Js ^s{n,T)u^ - - {J^^^^ + SJ^^J 

= J^' - ds It'' - s5u^ - ^ ( g„(5u^ + dg„w^ - dq, ^'^ 



In the above expression each of the terms inside the bracket vanishes because of thermodynamic 
identities (recall also that ^ = /is). 
It follows that 



A. 2 The Divergence of the fluid frame entropy current 



Using the relation that fJ-s = £, — y/^S.t^^.'^ the stress tensor current and the phase as given in ( pT^ ) 
can be rewritten as 

T^" = (pn + P^u^'u" + Pif + + 7^^" 

J'' - inu^ - + JLs 

We will find use, below, for the following thermodynamical relationships: 

_ Pn + P - mn 

T (A.5) 

dpn = p.dqn + Tds -~ Qsdfls 

Now divergence of su'' gives the following expression. 



^u.ds + s{d.u) 

= ^ [u.dpn + {pn + P) {d-u)] - j; [u.dqn + qn{d.u)] + j;U.dp.s 



(A.6) 



Using conservation of stress tensor (the equation to be used is u^d^T^^'^ = 0) one can evaluate the 



■^^ Recall that the symbol d denotes the change of a quantity under a frame change field redefinition, and that the 
microscopically defined field is taken to be the same in all frames so that, in particular, = 0. 
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expression [u.dpn + {pn + P) [d.u)]. 

U.dpn + {pn + P) {d.u) 



=fl [{u4){d.o + u-'{^.d)^ + {um-d) 



=^[p{d4) + u-'{^.d)^,] + p{Cd) 

Ms 



+ Pdiasd^ ( — 



(A.7) 



= ^ - (^) +Pdissd^ +U,d^'K^^ 

In the last line we have used the following identity derived from the curl-free condition. 



>u''{^.d)^, = ^{u.d){^^e) = -Ps{u.d)ps 



(A.8) 



Similarly using the current conservation equation one can evaluate the expression [u.9g„ + qn{d.u) 



u.dqn + qn{d.u) = ^{d.O + {Cd) 



d.Jdi 



(A.9) 



Adding all these equations one finds the following expression for the divergence of the entropy current 





'04- ps{^-d)ps 




T 



Ps_ _ (h 

pI p. 



(A.10) 



1 
T 



{d-Jdiss) 



The first line is zero if ps = PsQs and the second line can be written as 













.T. 



'^dis 



Pdi. 



9, 
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A. 3 Direct Thermodynamical determination of the modified phase entropy curent and 
its divergence 

As we have explained in the main text, a modified phase fluid description works in terms of a gradient 
vector defined by 

where = -pu'^ + = "M^t^ + (nf" 

In a modified phase frame, the stress tensor and charge current are taken to have the form 



T"^ = {ps + P)u''v!' + + /e^Co + 
where / = ^ = and all thermodynamical quantities are functions of T, /x, ^o- 

so ^0 



(A.I2) 
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(A.13) 



We will find use, below, for the following thermodynamical relationships: 

pn + P - 

T 

dpn = pdcin + Tds - f^od^o 

We now compute the divergence of the vector sm'* 
=u.ds + s{d.u) 

= ^ [U.dpn + {pn + P) {d.u)] - ^ [u.dQn + qn{d.u)] + ^U-dlls 

Using conservation of stress tensor (the equation to be used is Uvdfj,T^'^ = 0) one can evaluate the 
expression [u.dpn + {pn + P) {d-u)] as 



(A.14) 



(A.15) 



u.dpn + [Pn + P) {d.u) 

=^54/^0^) + fu,{ia.d)c^ + u,d^^^- 

Similarly using the current conservation equation one can evaluate the expression \u,.dqn + (7ri((?.u)] 

u.dqn + qnid.u) = d^if^^) - d^J^,^^ (A.16) 

The fact that the vector is curl free gives the following identity 

^o{u.d)^o + Un{^o-d)^o = -IJ-diss^Q{u.d)u^ + {C.d)iJ,diss (A.17) 

Using all these equations one finds the following expression for the divergence of the entropy 
current 



fj, ~ Uvd TT^" f f 



It follows that if we define an entropy current by 

then its divergence is given by 













T. 



T 



(A.18) 
(A.19) 
(A.20) 



In deriving this expression we have used the following equations. 

/ 



^Mu.d)^o + u^{^oMi;] 
f 



[-li'dissQ{.u.d)u^ + {C,-d)ndi8s] 



= - d„ 



rpC^f^di: 

da ( ICl^di. 



da 



~\~ f^diss 
H" f^diss 
H~ f^diss 
H~ l^dissP^ 



da 



da 



T 
T 
T 



j;^^{u-d)u. 
j;C{u.d)u, 



(A.21) 
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A. 3.1 Fluid frame entropy current transformed to the modified phase frame 



In this subsection we will check that the modified frame entropy current (A. 19) is precisely the current 
what we get by transforming the general fluid frame entropy current to modified phase variables. As 
we have explained in the main text, the relation between dissipative parameters in a general fluid 
frame and the corresponding modified phase frame is given by 



JLs = JLs + dqnu^ - die, - I t^d.ssW^ 



(A.22) 



where 'd' of any function denotes 



The transformation of the canonical fluid frame entropy current into the modified phase frame is 
given by 



diss 



T 



T 



dq„ 



T 



T 



"SO 



= s(eo)w^ 7f;Jd,ss - — d£,ou'' + y Co 



(A.23) 



s(eo)u^ 



T 



^ n jf" 

_ M .A- 



f f-^diss 
f f-^diss 



(a^^'^ + Co^) 



^ rj-^ diss rp 

In going from the second to third line we have used the following thermodynamic relation 

dpn = IJidqn + Tds - f^d£, 
In going from the third line to fourth line we have used the fact that 



f-^diss I ^ 



Co 



Using equation (A.23) and (A. 19) it follows that 



7P _ ta* 
jg — jg 



B. Details on Superfluid Thermodynamics from Gravity 

In this section we will use the bulk equations of motion to demonstrate that our solution obeys two 
Gibbs Duhem type relations. We will also review the demonstration of Q that the onshell action for 
the solution is the negative of its pressure, and that the equation for the infinitesimal variation of this 
pressure obeys a thermodynamical first law. 
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B.l Bulk equations, symmetries, conventions etc 

The Lagrangian we study is given by 



(B.l) 



)- 


r^3(#*)i 




2e2 



where V3{x) is taken to vanish at its minimum, = V — iAg, and Dg, = V + iAg, and all the potentials 

Vi, V2, V3 are real. 

The bulk equations that follow by extremizing this Lagrangian are given by 

Dg [F2(#*)5>*] = ^FabF'^" + ^Da<l>D-<l>* + ^ 

Rah - ^Rgab - 69ab = Tab = + Vb ^^''^^ 

where 

2Vl(#*) / e , i pec, p 

' ab — g2 \ ac-f 6T ^yab-^ ^ C1C2 

lab = ^ ( 2 -gabDcCpD (f) ) - ^ .^ gab 

Here small Latin letters denote the bulk coordinates ({r, .x, y, z}) and Greek letters denote the 
boundary coordinates ({i;,a;,y,2;}). Greek indices are lowered or raised using the metric r]. 

In this section we study solutions that preserve translational invariance in the field theory direc- 
tions. By scaling our solution we can always ensure that its horizon is located at r = 1; we make 
this choice in what follows. Our spacetimc has a d parameter sot of killing vectors that generate 
translations in the field theory directions. We now identify two special killing vectors among this set. 
Let fc" = {0,u''}, be the unique killing vector that is null on the horizon. Further let = {0,n'*}. 
Where n** is defined as 

pnvt: 

= iu = lim QaaA"- 

where A = A"- da is the gauge field and P^" = i]^^ + u^^u^ . 

Below we will identify the vector with the normal fluid four velocity, while the vector points in 
the spatial direction of the superfiuid velocity, in a frame in which the normal fluid is at rest. 
As we have explained, the most general ansatz for our solution is given by 

ds^ = —2g{r) Uf^dx'^dr — f{r) Uf^u^dx^dx^ + j{r) {u^^riv + u^n^) dx'^dx^ 

+ k(r)n^n^dx^'dx'' + r"^ P^.^dx^' dx" (B.3) 
A = A'{r)dr + H{r) u^'d^ + L{r) 

where 
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As we have mentioned above, we work in a gauge in which the scalar field is set to be real. This will 
be consistent with the equation of motion if the difference between the first and the second equation 
in (B.2) vanishes when (j)*{r) is set equal to One can check that when all potentials are real 

and there are no explicit dependence on the boundary coordinates the scalar field can be consistently 
chosen to be real if A'^{r) is set identically to zero. Consequently, on our solutions, 

A'' = and (f>{r) = (f>*{r) 

We now make the following convenient coordinate choices. Let v represent a flat normalized 
boundary coordinate in in the direction of constant vector u^^ and let x represent the coordinate is in 
the direction of constant vector n'^. In other words, in our new coordinate system, 

fc'^^jO, 1,0,0, •••} and r = {0,0, 1,0, •••} 

where we list a vector by (fc'', k"" , . . .). In this coordinate system the metric and the gauge field take 
the form 



ds'^ = 2g{r) dv dr - f{r)dv^ - 2j{r) dv dx + k{r) dx^ + dyfj 
A"- = 0, A" = iJ(r), = L{r) 



(B.4) 



Recall that in our coordinates the horizon occurs at r// = 1. Because the horizon is a null surface, 
the norm of the normal vector vanishes. As the oneform dual to the normal is simply dr we conclude 
that 

rj2(r)+/(r)fc(r) 



9"{r - 1) 



= (B.5) 

r=l 



g^{r)k{r) 

Recall also that, by definition, is a generator of the horizon, and so should be the dual to the normal 
dr at the horizon. This requires that gyy vanish at the horizon, so we conclude that /(r = 1) = 0. It 



follows then from ( [B.5| ) that j{r = 1) = 0. None of the other functions that appear in equation (B.4) 
are restricted at the horizon, except that they should be regular. 

At infinity we will impose the condition that our space is asymptotically AdS which implies that 



(B.6) 



and that 



lim g{r) 


= 1 + 


r^ 


lim /(r) 

r^oc 


= r'- 


f ^ + 0{r-') 


lim k{r) 

r— >-oo 


= r^- 




lim j{r) 

r—^OG 


- ii 


+ 0(r"4) 


lim Ay = 

r—>oo 


--^v = 


finite 


lim Ax = 




finite 



(B.7) 

lim Vi{<p(f>*)F'"'] = finite 

r— >oo 

lim [V5 Vi(#*)F^n = finite 



■^^The existence of such solution will impose some constraints on the potential such that the contribution of gauge 
field and matter field to the bulk stress tensor vanishes sufficiently rapidly as r oo. 
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The first two equations in (B.7) may be regarded as definitions of tlie constants and ^x- Tlie 
last two equations assert tlie finiteness of the boundary charge density and the boundary current in 
the X direction. 

B.2 Derivation of the Smarr Gibbs Duhem Relations 
B.2.1 Basic Strategy 

It is a geometrical fact (see e.g. section 5.3 of |5^) that if (^'^ is a killing vector 
It follows that, in our background, 

where we have used the Killing equation in the second line, the fact that V^C° is antisymmetric in the 
third line and the fact that all functions in our metric and gauge field depend only on r in the last 
line. Plugging in the explicit form of our metric and making different choices for the free indices we 
find 

= -ft [^,,'V'va-.l = -ft ( M<!i±™ ) 



V9RI = -dr [V55™5^'Vag,] = -9. 



r[k{r)fir)+f{r)] \ (B.9) 
9ir)VW) / 



„ r ^ ™ , n ^ ( r^[j{r)f'{r)- f{r)f{r)] 

= -Or [^/gg g Vah\ = -Cr - 



where = {0, 0, 0, 1, 0, • • • } 



The RHS of each of (B_^) are total derivatives. Our basic strategy is to find appropriate linear 
combinations of the four relations above so that the LHS of these equations also reduce to total 
derivatives onshell (i.e. upon using Einstein's equations). We will find three such linear combinations 
and so deduce the vanishing of three different expressions, each of which is a total derivative. We 
will then integrate these three equations to obtain three relations between quantities at infinity (i.e. 
conserved charges, currents, and superfluid velocities) and quantities at the horizon (entropy and 
temperature). Using these relationships we will be able to identify the vector (defined in terms of 
the generator of the horizon) with the normal velocity of the fluid, and also deduce two distinct Smarr 
Gibbs Duhem relations. In the next subsubsection we present the algebra involved in identifying the 
relevant total derivatives. In the subsequent subsection we will integrate these total derivatives to 
deduce physical conclusions. 
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B.2.2 Total derivatives that vanish onshell 



We use Einstein equations to simplify tlie LHS of tlie equations above. 



— Ih 



' n 



(B.IO) 



Here Tab is the bulk stress tensor. 

In order to do this we will now present and manipulate explicit expressions for the relevant bulk 
stress tensors. 

The contribution of electromagnetic and matter fields to the stress tensor above is given by 



-"afc — 'afc 

2Vi 



9ab 



CSC 'T-SC 

'-'ab — lab ^ 



^2(0) 



{daCl>db(l) + AaAbCl>'') + 



(B.ll) 



3e2 



gab 



where we have used the fact that the scalar field is real in the last equation. 

Using the fact that all components of the gauge field, metric and scalar field are functions of r 
only, one can simplify those components of the Maxwell equations where the free index is the boundary 
direction. 



*\pafj.] _ 



4 

V9 



(B.12) 



>dr.[VgVi{cj,)F^'^]^^A'^V2 



We now find explicit expressions (in terms of the functions that appear in our ansatz for the metric 



and gauge field) for the components of S^™ that we will need in the sequel. Using the equation (B.12) 



^(a,A,)F'^'- - ^{drA,)F- 



^A^dr iVgViWF^n - (V^i/i(0)F"-) 



(B.13) 



y^V,i(b) { -A,F 



3 



V5'^'V^2(0) 



Ia^a^ - ^ 

3 3 
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Similarly 



1 



(B.14) 



(3)e^ 



(B.15) 



V5 



3e2 



[A,A^ + A^A-] 



(B.16) 



In a similar manner we now find explicit expressions for the same components of 



V9 



V9V3W 
(3)e2 



3 

^3(0) 



(B.17) 



With these expressions in hand it is not difficult to determine three linear combinations of the equations 
appearing in ( B.lOj ) so that the RHS of of these equations is a total derivative. We find 
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[5" 



2 



(B.18) 



B.2.3 Integrating the Total Derivatives 



We now plug ( B.18 ) into (BJ) to obtain three total derivatives that vanish onshell. We now integrate 
these expressions from the horizon to infinity to obtain 



r[k{r)f{r) + j^r)] r^[k{r)f (r) + j{r)3' {r)] 



^[AAV9Vi{cl>)Fn]:ZT 



and 



r=l 



and 



4 [A {V9Vi(.<p)Fn]:zT 



The asymptotic expansion (B.4C) and (B.7) at infinity and the fact that 

/(r = l)=j(^ = l) = 



(B.19) 



(B.20) 



(B.21) 
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at the horizon aUow us to simphfy the LHSs of ( B.IE ), ( B.20 ) and ( B.21 ) as follows 



r[k{r)f(r)+f{r)] r^[k{r)f{r) +j{r)f{r)] 



9{r)^/W) 



-2/2 



23(1) 



r[k{r)f{r)+j\r)] r^[f{r)k'{r) 



g{r)^/k{r) 



2ko 



(B.22) 



r^W)nr)-f{r)nr)] 



= 2j2 



The RHSs of (|B.19 ), ( B.20 ) and ( B.21 ) may also be simplified upon noting that according to the 
prescription of AdS/CFT the boundary current J'^ is given by 



J'' = lim ( 4V5^i(<^)^^'' ) , = Hm g^^A'^ 



and also that 

Mr = 1) = - [f{r)H{r) + j(r)L(r)],.^, - 
(because / and j vanish at the horizon) and that 

F-'-(r - 1) = [V^A-],,^! = [/'■V,.A-],^i = 



because g'^^{r = 1) = Using these relations the RHSs of those three equations may be simplified as 
follows 



|[A„(v^ii(0)F-)];:r' 



(B.23) 



Finally plugging the simplifications ( |B.22D and ( [B.23D into ( |B.19D , ( |B.20D and ( |B.2lD obtain the 
final result of this subsubsection 



/2 
k2 

h 



45(1) 
167rG 



(B.24) 



4 

IBttG 



^vJx 
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B.2.4 Constraints on boundary stress tensor and current etc 

Let us first use the identities in (|B.24[ ) of tfie previous subsubsection to demonstrate tliat tlie boundary 
stress tensor dual to our gravitational solution takes the form 



T^u = {pn + P) U^Uy + Ps {Us)^{Us)u + PVp-v 



(B.25) 



In other words, u^, defined in this section as the killing vector that reduces to the black hole horizon, 
is also the normal fluid velocity, according to our definition of the normal fluid. In order to check 
the absence of a normal fluid super fluid cross term in (|B.25|) it is necessary and sufficient, in our 
coordinates, to check that 

' '^yy = ^ (B.26) 

Now using the usual definition of the stress tensor we find 

1 



T — 



T — 



-^yy ~ 



16nG 
1 

16nG 
1 

167rG 

T — 



-3/2 + 654 + 4fc2 
(-4j2) 



-/2 - 654 2 



(B.27) 



1 



16nG 



-h - 6.g4 - 4^2 - ^ 



Using these relations, the LHS of ( B.26 ) reduces to Using the last two equations in ( B.24 ) 
then immediately yields ( B.26 ). 

Let us now demonstrate the first relation in (|2.2|). In order to do this we note that the horizon 



area of our solution is equal to ■y/fc(l) and so that the entropy density of our solution is given by 
^'^'^^^ The periodicity of the Euclidean time circle in our circle is equal io P = ^ = An {jTjfj^ ■ It 



AG 

follows that 



Ts 



1 7fc(T)r(l) 



(B.28) 



167rG g{l) 

All other thermodynamical charges may be evaluated from the stress tensor as follows. In our 



coordinates 



Pn 



Ps — 

- P = T,, 



e 



T, 



yy 



T — 



(B.29) 



■^*In order to obtain these relations we worked, for concreteness, with a scalar field is dual to an operator of dimension 
2, i.e. one whose normalizable solutions look like 

4>2 



lim 0(r) = HI +0(r-*) 
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Substituting (B.27) in the first two equations of ( B.29 ) and tlien using ( B.24 ) and the last two 
equations of ( B.29 ) one finds tlie following relations 



167rG(P + p„) = — -y^^; h 167rGCi,q,i = -^^ — 777^; h leTrG^^^-u^jgn 



5(1) 



5(1) 



(B.30) 



(B.31) 



Finally using (B.28) these two equations turn into 

P + p„ =Ts + ^^M''q„ 

Ps = p-sqs 

the two Smarr-Gibbs-Duhem relations we set out to prove. 
B.3 The on-shell action 

In this subsection we will demonstrate that the onshell bulk action of our solution is the negative of 
its pressure. 

The bulk stress-tensor appearing in the Einstein equation that follows from the action in (5.1) is 
given by 

Tab = 72," + 7^6 + Qgab, 



(B.32) 



where 



' ah 



7'SC 
nh —~ 



V2{4>r) (Da4>Db(t>* + Da(^*Db(t> 1 



(B.33) 



2e2 



-9ab 



We consider solution of the form (5.5). Then it follows that the yy-componcnt of the stress-tensor 
is given by 



'^yy ~ 2 9yy ^) 



In deriving this relation we have crucially used the fact that the term Dy(j)Dy(j) in T^at and FycF'^y 
in Tmax is zero. This is because and the non-zero components of the gauge field is not a function 
of the coordinate y and the y-component of the gauge field is zero in our chosen form of the solution. 
Also we recall the identity for the Einstein Tensor E^jj, 



M 2 



Using the above identity and the Einstein equation Tyy = Eyy, we have from (B.34) 

9yy -J 

Now for the metric ansatz in (|5.5|) this expression reduces to 



(B.35) 



C 



1 



-g dr \r 



-9 9 



(B.36) 
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With the help of this relation our onshell action reduces to an integration of a total derivative 



^os — 



1 



IQttG 



d^x-r -V--9 .9' 



dr 



(B.37) 



Hence, after performing the above integration we get surface terms from the horizon and boundary. 
Since g^^ is zero at the horizon only the boundary term contributes and we have 



5, 



V0I4 



OS 



leTrG V r 



-V-5 9 



(B.38) 



bdy 



In order to avoid the singularities near the boundary we have to add the required counterterms to the 
on shell action 



CT 



IGttG 



bdy 



-7 (/C - 3 + 



(B.39) 



where 7 and K, are respectively the induced metric and extrinsic curvature of a constant r surface 
and is the mass counterterm for the scalar field. This mass counterterm is used to cancel any 
boundary divergence of the scalar field. In this subsection we keep it arbitrary as our argument here 
is independent of the specific value of v; we assume the scalar field to be norsmalizable. 

Let us now consider the following asymptotic expansion of the functions that appear in the metric 
and the scalar field 



lim g{r) 

r—^oo 

lim /(r) 
lim k{r) 
lim j{r) 

r-->oo 

lim 4>{r) 



= 1 



h 
h - 



J2_ 

^2 



+ 0(r-4) 



(B.40) 



With these expansions wc find that renormalized on shell action evaluates to 

1 



S 



ROS 



CT 



167rG 



(/2 + 6 54 + 4 fc2 - 2 A 



Now the boundary stress-tensor is given by 

1 



2n(bdy) 



lim Kj" - KSi", + 3(5'; 



5tt 



(B.41) 



(B.42) 



In the fluid dynamic limit the j/y-component of this boundary stress tensor yields the pressure of 
the fluid. Plugging in the metric ansatz ( ^.5| ) into ( |B.42 ) we have 

^ (/2 + 6 34 + 4 ^2 - 2 A 4) (B.43) 



7^(bdy) — p _ 

yy ~ i6^G 



Hence, comparing ( |B.43 ) and (B.41) we have 

Sros = -(V0I4) P. 

which implies that the renormalized on-shell action density is given by — P, as we set out to prove 

^^Here A is some constant fixed by the requirement that the limit in (B.42) is finite even when the non normalizable 
mode of the bulk scalar field is turned on. But when the scalar field is normalizable the limit exists for any constant A 
and in such cases the identity, proved here, is true irrespective of the value of A 
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C. Instabilities at large superfluid velocities 



In this section we study the linearized equations of superfluid dynamics in Herzog's model in a very 
simple context. First we restrict attention to the strict probe limit e — > oo. In this limit the only 
dynamical variables are and e. Next we restrict attention to perfect fluids. Finally we work at 
leading order in the e expansion. 

We assume that docj) = n{e) and dx4> = C where (j) is the superfluid phase. It is consistent, to linear 
order, to truncate to this special form (where we retain only two of the variables of fluid dynamics) if 
we also assume that e = e{x, t) and C, — t). More specifically we set 



e = eo + ^£6*'^''+''=^ 
C = Co + ,5Ce*"''+''=^ 

In order to obtain the dispersion relation for our fluctuations we need to solve the equations 



(C.l) 



9^J^ = 0, 9^e.-9.eA.=0 (C.2) 
where the second equation is nontrivial only when we choose (/.t, v) to be (0, 1). In the first equation 



is given by ( 6.17 ) retaining only terms of quadratic order in smallness of e and C, assuming that the 
two quantities scale in the same way. In other words the terms in the current that we have neglected 
are of the order ©(e^, C^, e^C; C^^)- 



The equations of motion have a solution of the form (C.l) only when the following 'characteristic' 
equation is obeyed 



epfc 
24 



6eofcCo+7eo" 3e5fc+24Co^ 
12 12 



= 0. 



(C.3) 



This equation yields the dispersion relation 



6 



Notice that uj develops an imaginary piece when ^ > \j 
we study has an instability when this inequality is satisfied. 

At the point of onset of instability the eigenvector that corresponds to the zero eigenvalue is given 

by 



Ir demonstrating that the system 



1 



1 1 

2 V 3 



(C.4) 



D. Limit C ^ 

As we have explained in the main text, the metric and gauge field dual to a super fluid flow must 



take the simplified form ( 7.21 ) in the limit C — >■ 0. In this Appendix we will explicitly verify that the 
results for our metric and gauge field at nonzero C reduce to this rotationally invariant form when C 



is set to zero, and read off all the unknown functions of radius in (7.21) 



Let us start with a metric and gauge field of the form (7.21) and rewrite it in the language of our 
metric and gauge field at nonzero That is we set C'' — Cn'^ where is a unit vector orthogonal 
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to u^, and eventually set C, to zero. The scalar vector and the tensor quantities defined above 
become 



2Si 



3 



2S'4 — Sq I Pi, 



-"6 / -> 



(D.l) 



Plugging these expressions into ( [7.2lD we find 



u^u^ + r P^^ 



dx^'dx'' 
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a^^dx^dx'^ 



Gi { {S4 + S6) + g2 ] Si 



Up^dx^dr 



Ti ( (^4 + 5'6)+ ) 5i 



UnUydx^dx^ 



^1 ( f ) ^3+V2 f ^ ) ^2 



Ufj^n^dx^dx'^ 



-5i 



(2^4 - ^6) 



P., 



dxf'dx" 



Vi ( [y4]. + V2 iy^y ] [Vi], 

ri(f) [F5]. + r2(f )[y2]. 



dxf^dx" 



dx'^dx'' + yT! f 7 ) [TiJ^.dx'^da;" 



(D.2) 



^1 - (^4 + 56) + ^2 - S 



dr 



_^ 1 
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+ 



^1 - (^4 + ^6) + H2 - Si 



£1 - 53 + £2 - ^2 
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We would now like to compare these expressions with the ^ — )• limit of the metric and gauge 
field at nonzero ^. There is, however, an immediate complication. The metric in (|D.2| ) is not in the 
gauge employed in our computation at nonzero C. In that more general calculation we have chosen a 
gauge such that the coefficient of P^^dx^dx^ term in the metric is r^. But in ( p.2| ) the coefficient of 
P^vdx^dx'' term is given by + where 



Kr) = 



(^) . , (^) 



^1 + 



-(254-56) 



So in order to compare (D.2) with the results at nonzero C we must first perform a coordinate 
redefinition upto first order in derivative. Define 



~2 2 

r — r 



such that 
r fa r ^ 

dr sa df 



(D.3) 



Mf)V 



Since this coordinate redefinition has a factor of it does not affect the gauge field as the gauge 
field is computed only upto order O (^)°- However the coordinate transformation does affect the 
metric, though only in the scalar sector (i.e. in the coefficients of scalar terms), basically because the 
coordinate transformation performed above is a scalar operation. Once this shift of r coordinate is 
implemented new terms are generated from the expansion of the uncharged black-brane metric (the 
metric at order C (^) )■ 

After performing the relevant coordinate transformation, the metric takes the following form 



ds-" 



-2g ( — I u'^dx^'dr 



2ru^ ( (u.d) (d-u) Ui,j+ VcF (^]au 



dxf'dx" 
dx^'dx" 

Mr)V 
2r 



u^dx^dr 



dxf^dx" 



Ufj^riiidx^dx^ 



(^) . , (^) 



-Si 



(254 - 56) 



h{f) 



rLaUudx'^ dx'^ 



(D.4) 



dx^dx" 



dx^'dx" + ^rii—] [Ti]f,^dx''dx'' S> + O 
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Now this metric matches the C 
among other conditions, that 

The constraints for the metric: 



hmit of the metric and gauge field in this paper provided, 



hm (5(72 {r) - 

hm (5/2 (r) = 
x-s-o 

hm 5K2(r) 

lim dJi(r) = 

x^o 

hm SJ^ir) - 
x^o 

lim (5J2(r) = 

x^o 

lim SKQ{r) 
x-s-o 



lim (5(/3(r) = lim (5(75 (r) — 

x->o x^O 

lim (5/3 (r) = lim (5/5 (r) = 

X-i-O x^o 

= hm Ji^3(r) = hm 6K5(r) = 
x^o x^o 

: lim 6J4(r) = lim dJr-,(r) = lim SJafr) = 

x^o x^o x^o 

: lim nW - ivi(r) 
x^o 2 



(D.5) 



lim Yi(r) 
x^o 



1 



lim 

x^o 



2 

SK^ir) 



lim SZi{r) 

X-i-O 



lim (5W5(r) 
x^o 



lim SKi(r) = 2 lim (5^2 (r) = r2(r) 

X^O x-i-O 

r/ie constraints for the gauge field 



lim SA^M 
x^o 


= lim (5yl6(r) 

x^o 




lim 5Hi{r) 
x-5-0 


= lim dHeir) 

x^o 


= Hi(r) 


lim (5i3(r) 
x^o 


= hm 6X4r) 

x^o 


= A(r) 


lim 5L2{r) 

x-^o 


= lim 6Xi{r) 





(D.6) 



It turns out all these constraints (and all others that are required to ensure that the bulk metric 
is an analytic function of the vector field w'' in the C ^ limit) are true upto the order the metric 
is calculated and using them one can read off Vi(r), V2(r), 7i(r) and T2{r) and hence the gauge 
transformation parameter h(r). Similarly the functions Giir), ^2(''), and J^2{f) appearing in 

Ufj^dx^^dr and Ufj^u^dx'^dx'^ terms are related to the corresponding functions appearing in the metric 
at nonzero C,. These relations also imply some constraint on the x ^ limit. But they are a bit 
complicated to write because they involve the function h(r) and its derivative with respect to r. Once 
Vi(r), V2ir), Ti{r) and T2{r) are all hxed, using them it becomes easier read off Gi{r), G2{r), J'lir) 
and J^2ir)- 
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